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Abstract 

From 1992 to 2007, HERA, the first electron-proton collider, operated at cms energies of about 
320 GeV and allowed the investigation of deep-inelastic and photoproduction processes at the 
highest energy scales accessed thus far. This review is an introduction to, and a summary of, the 
main results obtained at HERA during its operation. 
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1 Introduction 



HERA was the culmination of 50 years of experimentation with electron, later also muon and neutrino, 
beams to explore the structure of the proton, first in elastic and subsequently in inelastic scattering. The 
idea had been simple: use the pointlike lepton as a probe to study the strong interaction phenomena of 
the nucleons and their internal structure. HERA emerged from a series of electron-proton accelerator 
studies in the 70's as the most energetic version of an ep collider possible. It was built in the 80's for 
maximum possible luminosity, and with the capablility to scatter polarised electrons and positrons off 
protons, at a center of mass energy, ^/s, of about 320 GeV. The first HERA data were taken in summer 
1992. It ceased operations in June 2007 after a long, successful data taking period of 16 years. 

In deep inelastic lepton-nucleon scattering (DIS), the proton structure is probed by a virtual photon 
(7), or weak boson [Zq or W^), which carries a four momentum squared < s. The momentum 
transferred is inversely proportional to the spatial distance inside the proton which can be resolved 
with the photon or the weak boson probe. High energy DIS experiments are thus microscopes with a 
resolution power extending much beyond the eye, optical microscopes or modern laser techniques. The 
basic experimental interest and challenge over decades has naturally been to enlarge by increasing s. 

The ep collider at the DESY laboratory in Hamburg, Germany, had its first phase of operation, 
"HERA I", from 1992 through 2000. In this period, the collider experiments HI and ZEUS each 
recorded data corresponding to integrated luminosities of approximately 120 pb^^ of e~^p and 15pb~^ 
of e~p collisions. The energy of the electron (positron) beam was about 27.5 GeV. The energy of the 
proton beam was initially set to 820 GeV until it was increased to 920 GeV, in 1998, and kept at that 
value. 

The HERA collider was then upgraded to increase the specific luminosity by a factor of about four, 
as well as to provide longitudinally polarised lepton beams to the collider experiments. The second 
data taking phase, "HERA H" , began in 2003, after completion of the machine and detector upgrades 
and when unexpected large background problems were finally overcome. It ended in 2007. The HI 
and ZEUS experiments each recorded approximately 200 pb~^ of e~^p and 200 pb~^ of e~p data with 
electron (positron) energy of approximately 27.5 GeV and proton energy of 920 GeV. The lepton beams 
had an average polarisation of approximately ± 30% with roughly equal samples of opposite polarities 
recorded. 

In the last three months of HERA operation, data with lowered proton beam energies of 460 GeV 
and 575 GeV were taken, each experiment recording approximately 13pb~^ and 7pb~^, respectively. 
The primary purpose of this data is the measurement of the longitudinal proton structure function F^. 

The physics at an ep collider, which was to follow the fixed target experiments performed in the 
70's, was summarised in a well-known article by C. Llewellyn Smith and B. Wiik [IJ written 30 years ago. 
HERA has fulfilled those early expectations but delivered much more insight. While some predictions, 
such as the unification of electromagnetic and weak interactions, at values as large as the squared 
masses of the weak bosons, were confirmed, new, original and unforeseen, views on the structure of 
the proton and the dynamics of quark-gluon interactions were gained. The observations at HERA 
led to an enormous development in the understanding of the underlying field theory of the strong 
interactions. Quantum Chromodynamics (QCD). New fields of study emerged and existing fields entered 
a new precision phase; among these are the physics of high parton densities, the physics of high scale 
photoproduction and photon structure, the chromodynamic theory of heavy quark production and of 
hard diffraction. The limits set at HERA in the search for new heavy particles, as predicted, for example, 
in supersymmetric theories, are competitive with those set at the e~^e~ collider LEP and the pp collider 
Tevatron. As of this writing, no compelling indications of physics beyond the Standard Model have 
been found at these three colliders searching in the few hundred GeV energy range. The Large Hadron 
Collider (LHC) is now set to expand the search to the TeV scale. The proton structure determined 
with precision by HI and ZEUS at the ep collider HERA is of crucial importance [2] for the physics of 
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pp collisions at the LHC 

The current review is an overview of, and an introduction to, the physics results as obtained during 
the years of operation of HERA, with emphasis on the HERA I results since much of the more recent 
data is still being analysed. During the years 1992-2007, HI and ZEUS have each published more than 
one hundred and fifty journal publications. Given the space limitations, not all subjects can be covered 
in the present review. The reader is referred to the web pages of the experiments [3] and to the series of 
international workshops on "Deep Inelastic Scattering and QCD" [4J at which, since 1994, the physics 
at HERA and related results have been discussed annually, and in depth. This paper focuses on the 
results of ZEUS and HI and therefore cannot do justice to the impressive amount of theoretical work, 
on which the interpretations of the HERA data often rely. The readers are referred to the publications 
of the collaborations and papers cited therein. A review on early HERA results has been given in [5]. 
A recent general introduction to the physics of deep inelastic scattering is [6j. 

The paper is organised as follows: in Section 2 HERA accelerator issues are discussed, the ep 
kinematics and the two collider detector concepts are presented. Section 3 is devoted to the basic 
inclusive cross section measurements and the structure functions are derived. Section 4 presents the 
approaches of ZEUS and HI to extract the quark and gluon momentum distributions in the proton from 
QCD fits to the inclusive data. Section 5 discusses measurements of jet production and their importance 
to the determinations of the gluon distribution and the strong coupling constant, ag. Section 6 is devoted 
to the results of the production of charm and beauty quarks and their understanding in QCD. Section 7 
presents results on diffraction, including inclusive scattering, vector meson production and the new 
field of deeply virtual Compton scattering (DVCS). Section 8 summarises the results on tests of the 
electroweak theory, in neutral and in charged current ep interactions. Finally, in Section 9, the main 
results on searches for new physics at high scales are presented. 

The paper closes with a short summary and outlook. As this paper is being written there is still 
much analysis activity devoted to the publication of new and more accurate results, mainly based on the 
HERA II period of operation and upgraded detectors. In some areas, work on the combination of HI 
and ZEUS results has begun, which is expected to lead to the ultimate accuracy of HERA data. Some 
combined Hl-ZEUS data, such as those of inclusive DIS cross sections, have already been presented in 
preliminary form; however, as there are no combined published results available yet, their presentation 
is deferred essentially to a later review pubhcation. 

2 Accelerator and Detectors 
2.1 Introduction 

HERA, the "Hoch Energie Ring Anlage" at DESY, Hamburg, was the first, and so-far the only, accel- 
erator complex in which electrons and protons were collided. In the 60's and 70's several lepton-proton 
experiments were performed. Among these were the pioneering SLAC electron experiments [7J and the 
CERN muon experiments [HI [9] , in which a lepton beam of energy Ei between 20 and 280 GeV was 
scattered off a stationary nucleon target. The maximum energy squared in the center of mass system 
for fixed target lepton-proton scattering is s = 2MpEi with the proton mass Mp. The available cms 
energy sets the basic kinematic limits in any deep inelastic lepton-nucleon scattering (DIS) experiment. 

The SLAC experiment, in 1968, resolved the proton's structure down to a distance of 10~^^m at 
which it was discovered that the proton has point-like constituents. Since the energy available to lepton 
beams was limited to a few hundred GeV, several proposals were put forward to build an electron-proton 
collider as a two-ring structure with different types of particles being stored. In such a configuration 
the energy s becomes equal to AEgEp which is about 2Ep times larger than the fixed target energy, for a 
given lepton beam energy. Since the late 60's, Bjoern Wiik and colleagues had considered such machines 



5 



and proposed to probe proton's structure deeper with an ep collider at DORIS [TO], later at the e^e 
colliders PEP^, TRISTAN [12j and PETRA (PROPER) [13] and subsequently at the proton ring 
accelerator SPS at CERN (CHEEP) [T31|. Another ep collider proposal, to use the Tevatron accelerator, 
was put forward by a Canadian group (CHEER) [15]. When the decision to build a large electron- 
positron collider, LEP, at CERN was made, DESY endorsed the HERA proposal [16j. i.e. the plan to 
build an electron-proton collider at highest attainable energy, i.e. s ~ 10^ GeV^ with — 30 GeV and 
Ep ~ 800 GeV. This design was challenging since it was to use a chain of two times four pre-accelerators 
(see Table[l]) only some of which already existed at DESY, and enormous efforts were required to ensure 
its reliable operation and high luminosity. The design was made for positron- as well as electron-proton 
collisions and to achieve a large degree of longitudinal lepton beam polarisation, using the Sokolov- 
Ternov effect and spin-rotators around the interaction regions. 

Two general purpose collider detectors with nearly An acceptance were proposed in 1985, HI [18\ 
and ZEUS [19], and were built in the 80's. They were operated over the 16 years of HERA operation. 
Two further experiments at HERA were built and run in the fixed target mode. The HERMES exper- 
iment [20] (1994-2007) used the polarised e"^ beam to study spin effects in lepton-nucleon interactions 
using a polarised nuclear target. The HERA B experiment [21] (1998-2003) was designed to investigate 
B meson physics and nuclear effects in the interactions of the proton beam halo with a nuclear wire 
target. The results of the fixed target experiments are not discussed in this review. 

The collider detectors HI and ZEUS were designed primarily for deep inelastic ep scattering at 
highest momentum transfers, Q^, and large final state energies. Thus, much attention was paid to the 
electromagnetic and hadron calorimeters. The HI collaboration chose liquid argon as active material 
for their main calorimeter to maximize long term reliability. The ZEUS collaboration chose scintillator 
active media and uranium as the absorber material for the desired equalization of the calorimeter 
"e/vr" response to electrons and hadrons. The large calorimeters were complemented by large area wire 
chamber systems for muon momentum and hadron shower energy tail measurements. Because the e 
and p beam energies were very different, the detectors were asymmetric with extended coverage of the 
forward (proton beam) direction. Drift chambers inside the calorimeters, both in HI and in ZEUS, 
were segmented into a forward and a central part. Later, in HI starting in 1996 and in ZEUS from 2003 
onwards, silicon detectors near the beam pipe were installed for precision vertexing and tracking. Both 
apparatus were complemented with detector systems positioned near the beam axis in the accelerator 
tunnel, to measure backward photons and electrons, mainly for the determination of the interaction 
luminosity, and to tag leading protons and neutrons in the forward direction. Thanks to enormous 
efforts of several hundreds of engineers and physicists, both experiments took data for the entire time of 
HERA's operation at high efficiency. The only exception was a period of two years after the modification 
of the interaction region in 2001, which had caused unexpected problems to the ep operation. 

2.2 Accelerator 

Most of the data at HERA were taken with electron or positron energies of E^ ~ 27.6 GeV and Ep = 
920 GeV. In the first years, 1992-1997, Ep was set to 820 GeV. In the last months of operation, March- 
May (June) 2007, Ep was lowered to 460 (575) GeV. Typical beam currents were Ip ~ 100 mA and le 
ranged from about 35 mA at injection to 10 — 15 mA when the beam was dumped for refilling. Initially 
the use of PETRA as a pre-accelerator caused some problems due to eddy-current effects in HERA 
proton magnets at low energies and for tight injection aperture. With 174 colliding bunches the design 
luminosity of 1.4 ■ 10^^ cm~^s~^ was reached in 1997. The accelerator performance then stabilised and 
could be improved to give a total delivered integrated luminosity, between 1992 and mid 2000, of 
close to 200 pb~^ to each collider experiment. In the summer 2000, HERA operation was stopped for 
modifications of the interaction regions devoted to a further significant enhancement of the luminosity. 
For HERA, with the matched proton and electron beam sizes and head-on collisions, the luminosity 
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is essentially given by the product of the brightness, Np/ep, with the electron bunch current divided by 
the square root of the (3 functions in x and y. Methods considered to increase the luminosity included 
improvements in the injector chain which would overcome space charge effects limiting the brightness 
and an increase of currents with added power. Eventually, the luminosity increase was achieved by 
reducing the /3 functions at the interaction point (IP), from P^p = 7m and Pyp = 0.5 m to (3xp = 2.46m 
and f3yp = 0.28 m. This required the separation of the e and p beams much closer to the IP than in 
HERA I, i.e. at 11m instead of 24 m. Moreover, combined function (bending and focusing) magnets 
were introduced as close as 2 m to the HI and ZEUS interaction points. This forced HI and ZEUS to 
modify the inner detector configurations, and thus in 2000 the physics program for ep scattering with 
between about 0.1 and 2 GeV^ was terminated. It took some time until the effect of the superconducting 
magnets on the vacuum near the IR was finally understood and regular, approximately bi-monthly, 
warming up and cooling down cycles of these magnets were introduced which improved the vacuum. 

The period 2001-2003, following the modification of the IR, was extremely demanding. Beam induced 
background caused the central drift chambers of HI and ZEUS to trip at very low currents. Identification 
of the background sources and their suppression became of vital importance for the experiments and the 
luminosity upgrade to be successful. In an intense period of particularly strong collaboration between 
the machine and the detector experts, a series of beam based experiments and background simulations 
was performed over a period of more than one year [22l [23] . 

The different components of the background, those from the e and p beams, were investigated by 
varying the beam currents; special runs with one type of beam only were also made. Experimentally, 
it turned out that the dominant contribution was due to the dynamic interplay of synchroton induced 
effects with the proton beam: the new beam-line components combined with the bending of the electron 
beam closer to the IR than in HERA I caused ions to be released, which were scattered by the proton 
beam and its halo. Further problems arose due to the narrow aperture, from back-scattered synchrotron 
radiation and with leaks at flanges when the temperature load had risen too high due to mechanical 
or beam steering problems. The three successful counter measures were modifications of the beam 
absorber mask system, improved pumping and steady operation at slowly increasing currents to clean 
the surfaces using the synchrotron radiation itself. From the end of 2003 onwards, routine operation 
could be resumed with a specific luminosity increase as planned by a factor of about four and high 
positron/electron and proton currents. The annual delivered luminosity in 2006 was about 200 pb~^ , 
as much as had been delivered in the whole phase of HERA I. The machine was running very reliably 
and with high luminosity when its operation was ended. 

Table 1: Chain of (pre) accelerators at HERA and their energies. The main ring is 6223 m 
in circumference. High field (5 T) superconducting dipole magnets bend a proton ring of up 
to 1 TeV of energy. The electron ring requires high power, provided by a system of normal 
and superconducting 500 MHz cavities, to compensate for the synchrotron radiation losses. 
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Synchrotron emission leads to a build-up of transverse polarisation [17j with a characteristic time of 
about 30 minutes at HERA. Since 1994 this was used for the experimental programme of the HERMES 
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Figure 1: Luminosity delivered by HERA to each of the colhder experiments, HI and ZEUS, 
as a function of the number of days of operation. HERA was operated in an initial phase 
(HERA I), 1992-2000, with unpolarised lepton beams, mainly with . In the upgraded 
phase (HERA II), 2003-2007, the luminosity was increased and polarised data were taken 
with about equal amounts in terms of charge and polarisation states. The operation of 
HERA ended with two runs at 450 GeV (LER) and 575 GeV (MER) proton beam energy, in 
2007. 

experiment which was surrounded by spin rotators, a pair of dipole magnets flipping the spin from 
transverse to longitudinal orientation and back. In the first phase of HERA the collider experiments 
were not equipped with spin rotators. Polarisation effects, for unpolarised protons, result only from 
weak boson exchange at high as the interference of Zo-photon exchange in neutral current scattering 
occurs at at a size of order 10 ^Q^/GeV^. Thus, the spin rotators around HI and ZEUS were installed 
only for the luminosity upgrade phase, HERA II. For HERMES, depolarising effects from the coils of HI 
and ZEUS were compensated during HERA I with special magnets near the beam axis. Although the 
compensating magnets had to be removed when the focusing magnets were installed for the luminosity 
upgrade, HERA still achieved longitudinal beam polarisations with luminosity weighted means of 
typically 30% and maximum values of up to 40 — 50 %. 

Eventually both HI and ZEUS achieved total data collection efficiencies of 70-80% and collected 
large e^p data samples corresponding to integrated luminosities of nearly 500 pb~^ each, as is illustrated 
in Figure [H The HERA operations ended with efficient low proton beam energy runs, in which 13pb~^ 
at Ep = 460 GeV and 7pb~^ at an intermediate energy of 575 GeV were collected in only three months. 
This would have taken about a year in the old, HERA I, configuration. 

Physics at HERA would have profited from further running. Also one would have wished to accel- 
erate deuterons at HERA in order to study the structure of the neutron in the new kinematic range, as 
was proposed in 2003 [21]. Nevertheless, the 16 years of data taking and analysis so far have already had 
a great impact on the understanding of the partonic structure of the proton and on the development of 
the theory of strong interaction dynamics. After the final publications, expected in a few years' time, 
the HERA programme will have spanned nearly three decades of exciting physics and experimentation 
with the participation of perhaps a thousand experts. 
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2.3 Deep Inelastic Scattering Kinematics 

The deep inelastic ep scattering cross section, (Pa/dxdQ"^, of the inclusive reaction ep eX depends on 
the energy s = AE^.Ep and on two kinematic variables, usually taken to be the four-momentum transfer 
squared, Q^, and Bjorken x. At HERA, x is obtained from the measurement of the inelasticity and 
s as X = Q^/{sy). The salient feature of the HERA collider experiments is the possibility to determine 
the event kinematics in neutral current (NC) scattering from the electron e, or from the hadronic final 
state h, defined as all final state particles except the scattered lepton (and radiated photons associated 
with the lepton), or using a combination of the two. This leads to a maximum exploitation of the 
available kinematic range and to a redundant control of the measurement which was absent in the fixed 
target DIS measurements. Exploiting the NC calibration for the measurement of the charged current 
(CC) processes, in which there is no scattered electron, is an important means to reliably reconstruct the 
inclusive scattering kinematics for the reaction ep vX. The choice of the most appropriate kinematic 
reconstruction method for a given phase space region is a non-unique decision based on resolution, 
measurement accuracy and radiative correction effects. HI most frequently uses the "electron" and 
"sigma methods" while ZEUS most often uses the "PT" and "double angle methods" to reconstruct x 
and Q^. These methods are introduced below. 

The four-momentum transfer squared, Q^, can be calculated from the incoming (fc) and outgoing 
{k') electron four-momenta as follows 

Ql = -q^ = -{k - k'f = 2kk' = 2EeE[{l - cosae) = AE^E'^cos^ 6^/2 (1) 

neglecting the electron mass and defining 6e = tt — ae to be the angle between the scattered electron 
direction and the proton beam axis 0. The angle 9e is measured with the tracking detectors, or by 
the determination of the impact point at the face of the calorimeters and the event-wise reconstructed 
position of the interaction vertex. Here E'^ is the scattered electron energy, measured by the calorimeters. 
The inelasticity is given by k, k' and P, the four momentum of the incoming proton, and similarly to 
Qg, it can be calculated using the electron kinematics as 

Pa Pk' E' 

l/e = 7^ = l- ^ = l- TTTf (1 - cos^e). (2) 



Pk Pk 2E. 



e 



Often these relations are expressed as 



2 



^2 Pt,e 



in which Se is the difference between the energy, E'^, and the longitudinal momentum, p^ e, of the 
scattered electron, and pt^e is its transverse momentum. 

Similar relations are obtained from the hadronic final state reconstruction [25] 

^/i ^2 Pt,h /.X 



where 



is the hadronic E—pz variable and the sum, as for the transverse momentum, extends over the calorime- 
ter cells (and tracks in some cases) of the hadronic final state. In the reconstruction of the energy, the 



^ The angles are defined between the directions of the outgoing particles and the proton beam. This somewhat 
unfortunate choice can lead to some confusion since Rutherford 'back-scattering' at HERA is the scattering of the electron 
at high at small 'forward' angles. 
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momentum measurements from the tracking detectors are also used and the corresponding energy de- 
posits are removed from the sum. The hadronic energy is given by the relation 



Eh = (6) 
sm Oh 

The angle of the hadronic system, 6h, can be obtained from measurements of the energy and momentum 
components with the calorimeters as 

cos9h = ^ — ^'^ tan— = . (7) 

Plh + 2 pt^h 

In the naive quark parton model (QPM), 6h represents the angle of the struck quark. 

A straightforward calculation using Equations [T],[2] and Xe = Ql/si/e relates the uncertainties on the 
electron energy and angle to the uncertainties on and x as follows 

6QI/QI = 6E'jE'^ + tan{ej2)-6e, 

6Xe/Xe = 1/ye ■ SE'^/E'^ 

+ [(1 - ye)/ye ■ cot + tau {ej2)] ■ 56^. (8) 

The uncertainty of Bjorken x^ thus appears to diverge at small in the presence of a possible miscalibra- 
tion of the electron energy measurement or a misalignment affecting the polar angle measurement. On 
the contrary, the uncertainty of has no divergent behaviour, and indeed the electron method allows 

to be reconstructed optimally over the full kinematic range. There is also a term, oc tan (^e/2) ■ 56'e, 
which becomes sizable at very large angles, when the electron is scattered into the backward region, 
at small Q^. This was an important motivation for building dedicated Backward Detectors, EST (HI) 
and BPT (ZEUS), as supplements to the original apparatus, to extend the acceptance and to reduce 
the influence of the resulting cross section uncertainty to a tolerable level by measuring 9e extremely 
precisely at 6e approaching 180° (see Sections 2.4 and 3.4.2). 

Relations similar to Equation [H] can be obtained for the and x uncertainties in the hadronic final 
state reconstruction: 

5QI/QI = {2-yn)/{l-yH)-6Eu/Eu 

+ [2 cot eh + yh/{l- yh) ■ cot [9^/2)] ■ 66^ 
Sxh/xh = - yh) ■ SEh/Eh 

+ [-2 cot + (1 - 2yh)/{l - Vh) ■ cot {9^2)] ■ 59h. (9) 

Unlike the case of electron-based reconstruction, the uncertainties, both for and x, become large when 
y approaches 1. Thus the hadron reconstruction method and the electron method are complementary. 
It was proposed early in |26] to combine the variables Qg and yn in the "mixed method" reconstruction 
of the event kinematics in order to avoid the divergent feature of x^- 

A further refinement in the kinematics reconstruction was introduced with the "sigma method" |27j . 
Various versions of this method exist; the basic idea relies on the exploitation of the total E—p^ variable 

E-p, = E'^{1 - cos 9,) + Y.{Ei- p,,i) = Se + J^h, (10) 

i 

which for non-radiative events equals 2Ee. Thus yh, Equation HJ can be redefined and one obtains 

^2 Pt,e /„\ 

E -Pz ^-yT. 
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Scottered Electron Kinematics 



Hadronic Final State Kinematics 
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Figure 2: Electron scattering kinematics (left) and kinematics of the hadronic final state 
(right). Dashed: lines of constant scattering angle, de (left) and Qh (right), in the Q^, x plane. 
Dashed dotted: lines of constant energy, E'^ (left) and (right). For the electron energy 
a large region (left), in which E'^ ~ E^ is observed, because of the relation E'^ = Q"^ /AE^ + 
Ef.{l — y) ~ E'e, for small y and Q^. The E'^ distribution thus exhibits a characteristic 
"kinematic peak" at E^. which is often used to cross check or calibrate the scattered electron 
energy scale. One also recognises in the left plot that in forward direction, at very high Q^, 
the scattered electron carries very high energy of a few hundred GeV. The energy of the 
hadronic final state, E^, varies rapidly with Bjorken x at intermediate x values. At high x, 
the hadronic final state becomes very energetic and much of the energy is emitted in forward 
direction. At x = E^./ Ep = 0.03 the isoEh line is independent of and E^ = E^,. 



An important advantage of this construction is its reduced sensitivity to miscalibrations by using the 
final state energies and momenta instead of the incoming electron beam energy. The energy-momentum 
conservation requirement E—p^ — 2E^ has often been used to remove radiative events from the analysis. 

A further method, the "double angle method", developed independently in [28] and [29], uses the 
possibility to express and x by the electron and hadronic scattering angles which again have reduced 
sensitivity to energy miscalibrations 

tan(g,/2) 2 ^ 2 cot {ej2) 

tan (^e/2) + tan (^ J2) ' tan (ee/2) + tan (0 J2) ' ^ ' 

An application of this idea is the double angle energy calibration method, in which the electron energy 
scale is determined from the measurements of 9e and of 9^-, for regions of phase space in which both the 
electron and the hadrons are well contained in the apparatus. 

HI developed various versions of kinematic reconstruction methods. Diffractive analyses (see Sec- 
tion 7) often used a specific definition of the inelasticity as 

y = yl + yDA{i-yDA)- (13) 

A further idea has been to replace the electron beam energy in the calculation of x, similarly as it was 
replaced in This defines an extension of the E method which takes radiative effects at the lepton 
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vertex completely into account giving 



Xy = ^1 = (14) 

2Ep{E-p,)y^ 2EpJ:h ^ ' 

The hadronic final state, due to the particles lost in the beampipe as well as energy lost in the 
inactive material in the detector, is, in general, less well-determined than the scattered electron. In 
the "PT method" of reconstruction [30] used by the ZEUS collaboration, the fact that owing to the 
well- measured electron, 5pT = PT,h/PT,e gives a good event-by-event estimate of the hadronic energy 
loss is used to improve both the resolution and uncertainties of the reconstructed y and Q^. The PT 
method uses all measured variable, in a way specific to the ZEUS detector, to optimise the resolution 
over the entire kinematic range measured, namely. 

The variable 6pT is then substituted for 6^ in the formulae for the double angle method to determine 
y and Q^. The detector-specific function C is T,true,h/^h as a function of 9h, PT,h and 6pt and is 
determined from Monte Carlo simulations. 

The double differential cross section has to be independent of how and x are reconstructed. 
Many cross-checks of the results employ the underlying differences in the reconstruction algorithms and 
detector inputs. The measurement redundancies can also be expoited, as demonstrated in the case 
of the recent study of methods to average the HI and ZEUS measurements |2I] , to reduce the final 
measurement uncertainty to, eventually, about 1 %. 



2.4 Detectors 

The HI [22] and ZEUS [33] detectors are nearly Att hermetic apparatus with a solenoidal field of 1.2 T 
and 1.43 T, respectively, built to investigate high energy ep interactions at HERA. Their main design 
features are determined by the need to measure energies and momenta up to a few hundred GeV, as is 
shown in the kinematic plane plots (Figure [2]). 

Due to the large difference between Ep and Eg, the HERA apparatus, unlike similar devices at LEP 
or the Tevatron, are asymmetric. Large energy deposits, up to energies close to Ep, occur in the proton 
beam (forward) direction, and the particle multiplicity, due to the hard interaction and the proton 
remnant, is often very large. On the other hand, the hemisphere in the electron (backward) direction 
is less occupied, the energies deposited are below or comparable to E^ and the multiplicities are low. 
Such considerations have determined the final layouts of the HI and the ZEUS detectors, as illustrated 
with the event displays showing an NC event in HI, Figure [31 and a CC event in ZEUS, Figure jH 

The main component of the HI detector is the Liquid Argon Calorimeter (LAr). It has an inner 
electromagnetic part, which in the event of Figure [3] is seen to fully contain the scattered electron 
energy, and an outer hadronic part, which is seen to contain the fraction of jet energy leaking out of 
the electromagnetic LAr. The LAr is surrounded by a superconducting coil producing a solenoidal 
magnetic field of 1.2 T and an instrumented iron structure acting as a shower tail catcher and a muon 
detector. Tracks are measured in the central tracking detector which contains three silicon detectors, 
the forward, central and backward trackers, EST, CST and EST. These are followed towards outer radii 
by the Central Inner Proportional Chamber, CIP, and the Central Jet Chamber, CJC. The CJC has 
two concentric parts, separated by the Central Outer-z Chamber, COZ. Tracks in forward direction are 
measured in the Forward Tracking Detector, FTD, a set of planar drift chamber modules of different 
azimuthal orientation. 

The main component of the ZEUS detector is the Uranium-Scintillator calorimeter. Shown in 
Figure m are the forward, barrel and rear components (F/B/RCAL) which are surrounded by a muon 
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(jilt) Run 472542 Event 86273 RunPate 11/08/2006 




Figure 3: Deep inelastic event from neutral current ep scattering as registered in the HI 
apparatus in its upgraded, HERA II, configuration. The electron beam enters from the 
left and the proton beam from the right. On the left, an rz cross section shows the main 
components of the HI detector as denoted, see text. The event has a transverse momentum 
of the scattered electron of 91 GeV and a of 18600 GeV^ which correspond to y = 0.55 
and X = 0.33. Thus the electron scattering angle is about 9e = 30° and E'^ = 180 GeV, 
i.e. at high and large x the incoming electron is scattered in the forward direction and 
carries a rather large energy (cf. Figure[2]). The right top figure shows the event in the x,y 
projection exhibiting transverse momentum balance between the scattered electron and the 
hadronic final state as is characteristic for NC events. The lego plot visualizes the energy 
deposition in the LAr cells exhibiting the narrow jet structure of the hadronic final state 
emerging from the struck quark. 
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detector. The calorimeter is segmented into inner electromagnetic parts, of 5 x 20 cm^ tower sizes in the 
barrel, followed by hadronic parts of coarser granularity, 20 x 20 cm^. The space visible inside the BCAL 
is occupied by the coil giving a field of 1.43 T and a preshower detector. The tracking detector consists 
of a central cylindrical jet type drift chamber (CTD) and sets of planar drift chambers in forward 
(FTD) and rear (RTD) directions. During the upgrade phase, for HERA II, the FTD was modified by 
introducing straw tubes and inside the CTD a micro vertex detector (MVD) was added, comprising 
forward and central silicon strip detectors. 




ETA PHI 



UCAL transverse energy 



Figure 4: Deep inelastic event from ep charged current scattering as registered in the ZEUS 
apparatus, in its HERA I configuration. The electron beam enters from the left and the 
proton beam from the right. On the left, an rz cross section shows the main components 
of the ZEUS detector as denoted, see text. The right figure shows the measured transverse 
energy in the calorimeter. The event is a typical charged current scattering event with an 
energetic jet, unbalanced in transverse momentum, and some energy produced from the 
proton remnant in forward direction. For this event one finds 6h — 90° and pt^h — 33 GeV. 
From Equations HUT] this corresponds to ~ 2700 GeV^, y ^ 0.6 and thus x ~ 0.05, since 
at this time, HERA was operated with Ep = 820 GeV. 



Besides the main parts of the HI and ZEUS detectors, there were two instrumented areas further 
away from the interaction region, both in HI and ZEUS: 

• In the electron beam (backward) direction, several 'electron taggers' were placed, at 8 m, 33 m and 
103 m for HI and at 6 m, 8 m, 35 m and 44 m for ZEUS, in order to measure the scattered electron 
energy in a small range of de close to 180°. Some of these taggers were installed for only a part of 
the data taking period. The electron tagger system enabled, in different acceptance regions of i/e 
corresponding to the different positions, inelastic ep scattering processes to be tagged at very small 
Q^. The taggers thus allowed an efficient control and analyses of photoproduction processes. A 
process of particular interest is Bethe-Heitler scattering, ep — >■ ep^, whose cross-section is known 
theoretically to about 0.5% accuracy and is thus suitable for the determination of the interaction 
luminosity. Therefore both HI and ZEUS were also equipped with photon taggers positioned at 
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~ 100 m down the e beam line, at 180 degrees. The measurement accuracy of the luminosity thus 
far is about 1.5% for both experiments. 

• In the proton beam (forward) direction, protons and neutrons may be scattered at a small angle, 
via soft processes such as vacuum or pion exchange, respectively. Leading protons are a character- 
istic signature for diffractive processes such as hard diffraction or elastic vector meson production. 
Both HI and ZEUS installed so-called forward or leading-proton spectrometers for HERA I, EPS 
and LPS, at about 100 m distance from the interaction point. During the upgrade, HI kept its EPS 
and added another, very forward proton spectrometer, VEPS, positioned at 220 m away. ZEUS 
decided to remove the LPS during the upgrade phase. Similarly, HI operated a forward neutron 
calorimeter (ENC) throughout the HERA lifetime while ZEUS confined the ENC operation to the 
first phase of HERA. 

The detectors are complemented by systems of veto counters to reject beam background with time- 
of-flight measurements. In the forward direction, HI had a muon detector system comprising a toroid 
magnet inside a system of drift chambers. Similarly ZEUS had a toroid magnet with drift tubes for 
muon track reconstruction and limited-streamer tubes for triggering. 

The HI and ZEUS detectors functioned very reliably although intense and challenging efforts in 
the operation, calibration and maintenance were required. The main components, the central drift 
chambers, the big calorimeters and the large area muon detectors ran over 15 years with high efficiency 
and without a need for replacements of the basic hardware. Parts of the readout and trigger electronics 
were replaced or updated. In 1997 the HI central chamber had to be rewired when ageing effects were 
observed. Despite much higher backgrounds in HERA II, no ageing effects were seen in either of the 
HI and ZEUS tracking chambers, apart from a small reduction of pulse-heights with time. The purity 
of the argon liquid in HI was monitored and was found to worsen at a very small, tolerable rate of 
0.05 ppm per year. The argon was exchanged only once, in 1995. Similarly, the scintillator material in 
the ZEUS Uranium calorimeter survived all the operation well. 

The installation of focusing magnets near the beam pipe in the course of the luminosity upgrade, 
in 2001, required modifications of the HI and ZEUS interaction regions. It also provided the only 
significant time period during which large detector upgrades, as mentioned above, could be made. The 
HI detector upgrade program involved a large number of components; similarly ZEUS installed new 
hardware at various places of their apparatus. The original detector configurations are described in [52] 
and [33]. 



3 Proton Structure Functions 



3.1 Introduction 

The scattering amplitude for electron-proton scattering is a product of lepton and hadron currents times 
the propagator characteristic of the exchanged particle, a photon or Zq in neutral current scattering, 
a in charged current scattering. The inclusive scattering cross section therefore is given by the 
product of two tensors, 

5^-^i:.,if^T. (16) 

where j denotes the summation over 7, Zq exchange and their interference for NC, and j = or W~ 
for CC. The leptonic tensor L^" is related to the coupling of the electron with the exchanged boson 
and contains the electromagnetic or the weak couplings, such as the vector and axial- vector electron- 
couplings, Ve and Oe, in the NC case. This leptonic part of the cross section can be calculated exactly in 
the standard electroweak Ui x SU2 theory. The hadronic tensor, however, describing the interaction of 
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the exchanged boson with the proton, can only be reduced to a sum of structure functions, 
but not be fully calculated. Conservation laws reduce the number of basic structure functions in 
unpolarised ep scattering to i = 1 — 3. In perturbative QCD the structure functions are related to 
parton distributions / via coefficient functions C 

[^1,3, F,] =Y.t [1, ^]-C^l,2,3(-, ^, 4' «^(/^')) ■ /^(^' Z^?' Z^')' (17) 

where i sums the quark g, anti-quark q and gluon g contributions and fi{x) is the probability distribu- 
tion of the parton of type i to carry a fraction x of the proton's longitudinal momentum. The coefficient 
functions are exactly calculable but depend on the factorisation and renormalisation scales /i/ and /i^. 
The parton distributions are not calculable but have to be determined by experiment. Their de- 
pendence obeys evolution equations (see Section 4). A general factorisation theorem [31], however, has 
proven the parton distributions to be universal, i.e. to be independent of the type of hard scattering 
process. This makes deep inelastic lepton-nucleon scattering a most fundamental process: the parton 
distributions in the proton are measured best with a lepton probe and may be used to predict hard scat- 
tering cross sections at, for example, the LHC The parton distributions are derived from measurements 
of the structure functions in NC and CC scattering, as is discussed below. 

The following sections describe the structure function formalism, and present important NC and CC 
measurements. Both HI and ZEUS have collected data samples of roughly 50 — 100 pb^^ of unpolarised 
e^p data and of similar luminosity in the two charge and opposite helicity states (see Table[2]). 



Table 2: Integrated luminosity (L), in pb~^, collected by HI and ZEUS, for specific data se- 
lections, distinguished according to the beam charge and helicity (P) state. The unpolarised 
data were taken in HERA I. The Hl-ZEUS differences in polarisation and luminosity are due 
to somewhat different analysis criteria, data taking efficiencies and an asymmetry in rotator 
switches between HI and ZEUS during the e~p run, in the spring of 2006. While this review 
is being written a significant part of the polarised data, taken in the years 2004-2007, is still 
being analysed. 
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3.2 Structure Functions and Parton Distributions 

The neutral current deep inelastic ep scattering cross section, at tree level, is given by a sum of gener- 
alised structure functions according to 

where the electromagnetic coupling, the photon propagator and a helicity factor are absorbed in the 
definition of a reduced cross section cr^, and Kj. = 1 ± (1 — y)^. The functions F2 and XF3 depend 
on the lepton beam charge and polarisation (P) and on the electroweak parameters, e.g. on the weak 
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boson masses {Mz, Mw) and the fine-structure constant a, as [35]; 

F± = F2 + ^zi-Ve T Pae) ■ Ff + kI{vI + al± 2Pvea,) ■ 
xFf = Kz{±ae + Pve)-xFf + K\{T'2veae-P{yl + al))-xFi, (19) 

with K2(Q2) = Q'^/[{Q^ + Ml){4:Sin^ 6 cos^ 6)] and the weak mixing angle sin^ 6 = 1-M^/M|. In the 
hadronic tensor decomposition [SB] the structure functions are well defined quantities, independently to 
the parton model. In the Quark Parton Model (QPM) the longitudinal structure function is zero [3Z] 
and the two other functions are given by the sums and differences of quark (g) and anti-quark (g) 
distributions as 

{F^,Ff,F^) = xY.{el2e,v,,vl + al){q + q) 

{xFf.xFi) = 2xJ2ie,a,,v,a^)iq~q), (20) 

where the sum extends over all up and down type quarks and = e„, denotes the electric charge of 
up- or down-type quarks. The vector and axial- vector weak couplings of the fermions (/ = e, u, d) to 
the Zq boson in the standard electroweak model are given by 

f/ = z/ — e/2 sin^ 6 af = if (21) 

where e/ = —1,2/3,-1/3 and if = I{f)z,L = —1/2)1/2,-1/2 denotes the left-handed weak isospin 
charges, respectively. Thus the vector coupling of the electron, for example, is very small, Ve = —1/2 + 
2 sin^ B ~ 0, since the weak mixing angle is roughly equal to 1/4. In a famous experiment [38] in 
polarised electron-nucleon scattering at very low Q^, it was observed that a^. ~ —1/2 which proved the 
absence of right-handed weak currents in the available energy range, and thus the validity of the most 
simple weak doublet structure, now called the standard electroweak theory [39] . 

Out of the eight structure functions entering in Equation [18] four have been accessed at HERA. 
At low Q^, the cross section is determined by F2 and F^. The longitudinal structure function Fl is 
measured at low x. The large gluon density in this region leads to a rather large value of F^ compared 
to those from the fixed target measurements which were measured at a much larger x. 

At low and low y the reduced cross section to a very good approximation is given by cxr = 
F2{x, Q^). Indeed the accurate knowledge of the proton structure, at low x from HERA and at large x 
from fixed target experiments, is mostly due to precision measurements of F2. At y > 0.5, F^ makes 
a sizeable contribution to Cr^NC- In the DGLAP approximation of perturbative QCD, to lowest order, 
the longitudinal structure function is given by 
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(22) 



which at low x is dominated by the gluon contribution (see below). A measurement of Fl requires a 
variation of the beam energy since ar,Nc{x, s) = F2{x, Q"^) — y'^Fi^x, Q'^)/Y^; Fl is obtained as the 
slope in y'^/Y^ at fixed x and [y = Q'^/sx). Two further structure functions can be accessed with 
cross section asymmetry measurements in which the charge and/or the polarisation of the lepton beam 
are varied. 

Defining a reduced cross section, similar to ar,NC, for the inclusive polarised charged current e^p 
scattering as 



2vrx 

" ' ' dxdQ2 

one obtains a sum of structure functions, analogous to Equation [TH as 



dVcc 



(23) 



<cc = ^(^2^ T ^xWt - ^Wt) . (24) 



17 



Here the Fermi constant is defined [JT] using the weak boson masses, in very good agreement with 
Gp determined from the measurement of the muon hfetime 02]. In the QPM (where = 0), the 
structure functions represent sums and differences, that depend on the charge of the lepton beam, of 
quark and anti-quark distributions given by 

= x(U + D) , xW^ = x{D - U) , = x{U + D) , xW.^ = x{U - D) . (25) 

Using these equations one finds 

a+cc-xU+il-y^xD, (26) 

a-cc-^U+il-yfxD. (27) 

Combined with Equation [191 which approximately reduces to 

<ivc - [cu{U + U) + Cd{D + D)] + Kz[du{U -U) + d,{D - D)] 

with Cu,d = el ,^ + Kz{-Ve =F Pae)eu,dVu,d and du,d = ±aeau,deu,d, (28) 

one finds that the NC and CC cross section measurements at HERA determine the complete set U, D, 
U and D, i.e. the sum of up-type, of down-type and of their anti-quark-type distributions. Below the 
b quark mass threshold, these are related to the individual quark distributions as follows 

U = u + c U = u + c D = d + s D = d + s. (29) 

Assuming symmetry between sea quarks and anti-quarks, the valence quark distributions result from 

u^ = U-U dy = D-D. (30) 

Equation [28] also shows that the vector and axial vector couplings of up and down type quarks can be 
determined at HERA. As will be discussed below, recent QCD fits have given quite remarkable results 
on joint determinations of parton distributions and of the weak neutral current coupling constants of 
the up and down type quarks. One notices that the underlying assumption here is that all up as all 
down type quarks have the same couplings, e.g. a„ = ttc, which holds in the Standard Model. 

3.3 Measurement Techniques 

The inclusive neutral current scattering process is identified by the unique signature of an isolated 
electron cluster above some minimum energy, of a few GeV, in the LAr/Spacal (HI) or the U/Sc 
(ZEUS) calorimeters. In charged current scattering, a large transverse momentum, summed over all 
hadronic final state particles, is required (e.g. prih) > 12 GeV). The NC and CC scattering events 
are required to have an interaction vertex within the nominal interaction region, ±35 cm using the 
central trackers or, at the trigger level, of ±2ns from the nominal interaction time. Unlike HI, the 
time resolution of the ZEUS calorimeter is suitable for distinguishing genuine ep scattering from beam 
background outside the interaction region. 

NC events, unlike CC events, are balanced in the transverse plane. Therefore the hadronic calorime- 
ter can be calibrated requiring the ratio Pr(^)/PT(e) to be close to one. The calibration of the elec- 
tromagnetic calorimeters is based on several types of events, in which the electron energy scale can 
be established, such as the QED Compton process, elastic p production in DIS and events within the 
"kinematic peak" region of phase space, in which the scattered electron has an energy very close to 
the initial beam energy. Also, comparisons with the tracking detector measurements and the electron 
energy determined with the double angle method (see Section [273l) are used. Further adjustments to 
the linearity of the energy scale determination are made, for example in the case of the HI Spacal 
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calorimeter, with the reconstruction of the ttq and the J /if) mass. The calorimeters are aligned by 
matching the energy deposits with tracks found by the tracking system; also elastic Compton events, 
with characteristic back-to-back e — 7 energy deposits in the transverse plane, are exploited for this 
purpose. 

The basic measurement technique is to compare the reconstructed event distributions with complete 
Monte Carlo simulation samples with much larger statistics than the data. Backgrounds are subtracted 
from the data; most important of these are the events due to photoproduction where the scattered 
electron escapes near the beam axis and an energy deposit in the calorimeter is wrongly attributed to 
be the scattered electron. At small scattered electron energies, E'^ > 2 GeV, this background may well 
exceed the genuine DIS contribution. Therefore, in some analyses the calorimetric electron identification 
is complemented by a track associated with the calorimeter cluster. Additionally, efficiencies, such as 
those from the trigger condition used to select the CC or NC events, are determined from the data and 
are used to correct the simulation. Radiative corrections to the tree level cross section [H] in NC are 
much reduced by requiring the longitudinal momentum in the event to be conserved, i.e. E — pz — 2Ee 
(see Section 12.31) . The simulation is normalised to the measured luminosity, which is known to within 
1 — 2 % at HERA using Bethe-Heitler scattering, ep epy. The cross section, binned according 
to detector resolution, is then determined as the ratio of the background subtracted number of DIS 
events in data to the simulated events, each corrected for the luminosity of the data sample and the 
simulation, respectively, and multiphed by the DIS cross section used in the simulation. With such 
techniques, refined in several aspects, both ZEUS [13] and HI [H] have published, since 1993, about 15 
papers each on the measurements of NC and CC cross sections at HERA. 

3.4 Low and x Results 

3.4.1 The Discovery of the Rise of F2(x, Q^) 

The first measurements of the proton structure function E2{x,Q'^) , shown in Figure El were based 
on only 0.03 pb~^ of data, taken in 1992. With these first measurements the rise of F2 towards low 
X was discovered. This rise is in agreement with general expectations on the low x (large 00 = 1/x) 
asymptotic limit of QCD [45j; however, QCD could not predict the actual scale (Q^) at which the limit 
would be applicable. The fixed target data available up to 1992 led to many extrapolations towards 
low X] for example, classic Regge theory ^46j expected F2 to continue to be fiat at low x. The curve 
shown in Figure was the expectation of a model [17], in which the parton distributions were generated 
dynamically and evolved from a very low scale ~ 0.3 GeV^ towards higher Q^. In this approach a 
rise towards low x was indeed expected. 

With increasing luminosity and improving understanding of such measurements, the accuracy in- 
creased from about 20% in the first data to 2% for the most recent measurements. The data of HI 
and ZEUS have been rather consistent and mostly in good agreement for the various NC and CC 
cross section data. The best published measurements, thus far, of the reduced cross section at low 

(< 200 GeV^) are displayed in Figure [61 The left plot demonstrates the rise of the cross section 
from the rise of F2 towards low x in all bins of Q^. In the intermediate range the cross section, 
ar = F2 — y'^Fi/Y^, is observed to fiatten or even turning over under the infiuence of the Fl term at 
high y. Dedicated measurements and analyses are being performed to extract F^ as is discussed below. 

Another dramatic observation on the behaviour of F2{x,Q'^) was its strong dependence at fixed x 
on Q^, (see Figure [6] right). Thus the very first observations of scaling, the independence of F2{x,Q^) 
on Q^, has been confirmed only for x ~ 0.1 where it was first discovered at SLACB The HI and ZEUS 

^ For rather accidental reasons the energy and angular range of the SLAC ep spectrometer was such that the mea- 
surements correspond to x ~ 0.1. Therefore the first observations of inelastic DIS [7 occurred in a region where the 
behaviour of F2 could rather directly be related to the quark interpretation of the proton structure. Departures from this 
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Figure 5: The first measurements of HI (solid points) and ZEUS (open points) of the proton 
structure function F2{x, Q"^) based on the data taken in 1992 shown as a function of Bjorken 
X. The HI data at low (larger) x were obtained with the electron (mixed) method. The 
ZEUS data were obtained with the double angle method (see text). The HERA collider is 
equivalent to a fixed target lepton proton beam experiment of about 50TeV lepton beam 
energy. The HERA experiments were thus able to extend the kinematic range of the F2 data 
provided by the fixed target electron (SLAG) and muon (BCDMS, NMC) proton experiments 
by two orders of magnitude into the then-unknown domain of low x. The curve shows a 
result of the so-called dynamic parton distribution approach. 

measurements illustrate the important differences between the scaling violations at high x, where F2 
decreases as a function of due to gluon bremsstrahlung, and at low x, where F2 increases as a 
function of mainly due to qq pair production from photon-gluon fusion. The strong rise with at 
low X is a major discovery of HERA. 

3.4.2 Remarks on Low x Physics 

Confrontation of the measurements with the DGLAP evolution equations (see Section 4) leads to the 
result that at low x the quark contribution to the InQ^ derivatives is nearly negligible. This has two 
consequences: the measurement as presented in Figure [6] (right) can be interpreted as the observation 
of a very large gluon distribution xg{x,Q'^) at low x which at fixed rises with x. Furthermore, 
neglecting the quark part (see [H]), one can infer from the second DGLAP equation that xg{x,Q^) oc 
exp^c InT ln(l/a;) with T = <ys{Q^)/cis{Qo) [SQ] and get a valid estimate of F2 and its behaviour 
in and x. Indeed one is thus able to obtain that F2{x,Q'^) rises towards low x as a;~^ with A oc 
\J\tiT / ln(l/a;). The rise is thus expected to be stronger as increases as indeed is observed (see 
Figure [6] left). A numerical study of this behaviour was performed in [51]. Much deeper theoretical 
evaluations can be found in [52]. Apart from the numerical success of this and further studies (see for 
example [53]) one finds that the proton structure at low x differs qualitatively from the large x limit: 
while the latter is given by the valence quark behaviour, the low x region is dominated by the dynamics 
of the QCD vacuum, the gluons. 

behaviour, which were expected in field theory, were first observed at Fermilab [15], and with more accurate data, they 
could be shown subsequently to be logarithmic in . This was a major success in the process of establishing Quantum 
Chromodynamics as the appropriate field theory for quark-gluon dynamics. 
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Figure 6: Deep inelastic neutral current e~^p scattering cross section measurements, from the 
HERA I data taking period, as a function of x in different intervals (left) and for three 
selected x bins as a function of (right). Open points: HI, open squares: ZEUS. The error 
bars show the total uncertainty. The curves are NLO QCD fits as performed by HI and 
ZEUS to their own data. The flattening of the curves for the two lower x values at large 
is due to the effect of Fl (right plot). At x = 0.25 the data and the curves decrease visibly 
for above a few 1000 GeV^. This effect results from the 7Z interference which reduces 
the cross section in positron-proton scattering at high Q^, cf. Equation [19] , oc a^xF^^ with 
ae = -1/2. 



Studies quoted above, and QCD fits discussed below, appear to validate the application of DGLAP 
equations over most of the HERA kinematic range. However, questions remain; because x is as small 
as 10~^ for > IGeV^, there are expectations that terms neglected in DGLAP evolution, those 
proportional to ln(l/x), should become sizable and important [54j- Also, as becomes small, the 
applicability of perturbative QCD becomes suspect and impossible when as approaches 1. 

The rise of F2 towards low x has been examined quantitatively [55l[56] (see Figure [7j). For x < 0.01, 
where the valence quarks have a negligible role, the rise at fixed may be expressed as F2 = Cx^^. 
It has been measured that for > 3.5 GeV^, i.e. in the DIS region, the parameter C is approximately 
constant and A = 0.048 ■ (InQV^^) with A = 0.29 GeV. At lower the rise becomes weaker and C 
dependent on Q^. This may be interpreted as a change from the partonic behaviour, at larger Q^, to a 
hadronic behaviour [57] at lower Q^, the transition taking place near ~ 1 GeV^ corresponding to a 
resolution distance of about 0.3 fm. The so-called BPT data [58] below 1 GeV^ which were obtained by 
ZEUS with a dedicated detector placed near to the beampipe, together with the precision measurements 
based on the central apparatus of HI and ZEUS, play a crucial role in this observation. 

The transition region is often studied in terms of the "colour dipole model" (CDM), which will also 
be discussed in connection with diffractive DIS in Section 7. The CDM is able to provide predictions, 
using a phenomenological ansatz for a hard scattering cross section, for a variety of for the 

structure functions F2, Fl but as well for the diffractive and heavy quark structure functions, which 
are introduced below. The strict relations between these predictions lead to crucial tests of the validity 
of this approach [60] , results of which can still be expected at HERA. The CDM model fundamentally 
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Figure 7: Rise of F2 versus x as a function of Q^, see text. Note that with Q"^ also x varies, 
roughly as x ~ Q^/O.ls. The figure is from 



does not distinguish between the partonic or hadronic behaviour of the proton and thus is applicable 
for large as well as very small values of in DIS. It is to be also noted that the transverse size of the 
dipole is 2/ a/Q^; thus by varying the momentum transfer, one is able to probe the transverse extension 
of the proton. Most notably at HERA, the transverse size of gluons was measured for the first time 

The HERA physics of low x has important relations to other fields. The knowledge of the rising 
parton densities is of importance for investigations of the quark-gluon plasma phase in nucleus-nucleus 
collisions, as are studied at RHIC [62] and at the forthcoming LHC, as well as for cosmic super-high- 
energy neutrino physics [63]. Furthermore there is a wide range of forward scattering at the LHC, in 
which low X values are involved and it is, so far, not clear whether one may use the DGLAP equations, 
evolved to the LHC range, in order to understand large rapidity particle production characteristics. It 
is possible that modifications from the large contribution of ln(l/x) terms in this domain are needed 
for a full description. 



3.4.3 The Longitudinal Structure Function 

There is a strong interest in Fl since the longitudinal structure function is a measure of the gluon distri- 
bution which is independent and complementary to the derivative dF2{x,Q'^) / dhiQ'^ . A measurement 
of Q^) represents a crucial test of QCD to high orders as had long been studied experimentally fBl] 
and theoretically ^3]. The importance of this test at HERA rests also on the observation that the de- 
terminations of the gluon density at low x and low (< 10 GeV^) are not constrained well because of 
the lack of range and the uncertainties of the measurements of dF2{x,Q'^) / dhiQ^ . 

At lowest X, corresponding to large y, the reduced cross section departs from F2 because of Fl- This 
had been seen for the first time at HERA in [66]. A subsequent accurate measurement of the cross 
section up to high values of y is shown in Figure [HI This measurement is particularly difficult because 
at high y — 1 — E'^/E^, the scattered electron energy is rather small; E'^ = 2.7 GeV for y = 0.9. At such 
low energies the DIS electron signal is mimicked by hadronic energy depositions in photoproduction 
processes which poses severe difficulties in extracting the genuine DIS contribution. 

The measurements of the high y cross section have been used by HI to determine the longitudinal 
structure function at low x. This is enabled by the large kinematic coverage of HERA: the range in 
y extends over two orders of magnitude, thus at lower y (< 0.35) one is able to determine F2 which 
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Figure 8: Left: Measured reduced cross section as a function of x in the mid region. The 
data at lowest x are seen to depart from the extrapolated behaviour of F2. This departure, 
for each bin, occurs at about y ^ 0.5 and is attributed to the effect of the longitudinal 
structure function Fl which at low x is large. The curves are the NLO QCD fit result on 
for Fl = (dashed) and for Fl (solid) obtained in the fit. Right: The longitudinal structure 
function Fl{x, Q^) for different bins as obtained by HI at low x, and by charged lepton- 
proton scattering experiments at large x (see [67]). The curves are from the HI NLO 
QCD fit, described below, and its uncertainties. 

then may be extrapolated to the high y region. At larger Q^, in the DIS region, as is illustrated in 
Figure [HI one may use the NLO QCD fit to predict F2 in the region where only (Jr is actually measured 
and subtract the F2 part from the cross section. At lower Q"^ pQCD cannot be trusted and different 
means have been developed by HI [671 EH] to simultaneously determine Fl with F2. The published HI 
data on Fl are shown in Figure [HI further results were presented subsequently in [6H] • The principal 
problem of such determinations is of course their unavoidable dependence on the knowledge of F2. 

During the last 4 months of its operation, HERA was operated at reduced proton beam energies. 
Data samples of 13pb~^ at 460 GeV and of 7pb^^ at 565 GeV were collected by both ZEUS and HI. 
This data will allow a model-independent separation of the contributions of F2 and Fl- Since the 
luminosity behaves as L oc -E^, a reduction to half the energy requires 4 times longer running. A more 
precise measurement for statistical reasons will thus be restricted to the lower region. As this paper 
is being written, both the ZEUS and HI Collaborations are analysing the low energy run data and 
results are much anticipated. These will complement and perhaps validate previous attempts by the 
HI Collaboration to extract Fl as were summarised briefly above. An attempt will also be devoted to 
the extraction of the diffractive longitudinal structure function, F^ , which is of particular theoretical 
interest and about which nothing is known experimentally so far [HS] • 

3.5 High Q2 Results 

The large cms energy of HERA has permitted the extension of the measurements of DIS lepton-proton 
scattering cross sections, as compared to previous fixed target experiments by about two orders of 
magnitude towards high momentum transfers squared, Q^. This can be seen in Figure [9] (left) which 
presents the cross section measurements as obtained previously in fixed target experiments and those 
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from HI and ZEUS as a function of for the full range of Bjorken x covered. One observes scaling 
to hold for the complete range at x ~ 0.1. At larger x, negative scaling violations, due to quark 
bremsstrahlung, continue to exist up to high Q^. One also observes in this range the tendency of 
the cross section to decrease stronger at highest which marks the onset of the effects of the 7Z 
interference contribution which for e'^p scattering is destructive, (see Equation [19]) . Figure [9] (right) 
shows the HI and ZEUS data in the high region as a function of x. The data agree rather well and 
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Figure 9: Deep inelastic neutral current e^p scattering cross section data from the HERA 
I data taking period. Left: NC reduced cross section data versus for different x. The 
measurements are consistent with earlier DIS fixed target measurements but considerably 
extend the kinematic range in and x. Right: NC cross section data versus x in different 
intervals of . The curves are NLO QCD fits as performed by HI and ZEUS to their own 
data (see Sections 4 and 5.4). 

are described by NLO QCD fits. HERA has thus established perturbative QCD to be valid up to 
values beyond 10^ GeV^ which was one of the questions posed initially pLj when one envisaged a high 
energy ep collider to be built. 

Unlike fixed target experiments, which required separate electron/muon- and neutrino- nucleon scat- 
tering apparatus, neutral and charged current data could be taken simultaneously at HERA because 
at larger > 100 GeV^ the weak interaction process ep vX becomes measurable. Figure [TUl shows 
the published results from ZEUS and HI on e^p CC scattering. These data are of unique nature: they 
further support QCD, allow flavour contributions to proton's structure to be disentangled (Equation [25!) 
and, as is discussed below, they demonstrate clearly that the electromagnetic and the weak interactions 
become of similar strength for ~ ^wz- 

The data presented here belong to the initial (HERA I) period of data taking. The luminosity 
increase and the polarisation in HERA II have much enhanced the importance of the high NC and 
CC data, which will be published in the future. Yet, the data from the HERA I period have already 
had a major impact on the determinations of the quark and the gluon distributions in the proton, and 
also on the determination of the strong coupling constant as is discussed in the following. 
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Figure 10: Deep inelastic charged current e^p scattering cross section for between 300 
and 30000 GeV^, measured by HI and ZEUS using data taken in HERA I. The errors shown 
are the total errors, which are dominated by the statistical uncertainty. The curves are NLO 
QCD fits as performed by HI and ZEUS to their own data (see Sections 4 and 5.4). 



4 QCD Fits 



4.1 Introduction 

In perturbative QCD the quark and gluon distributions D{x,Qq) = q,g, determined at an initial value 
Ql, can be evolved further in using a system of coupled integro-differential equations (DGLAP 
equations) 
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which was introduced independently by Dokshitzer |70j, Gribov and Lipatov [7T] and by Altarelli and 
Parisi ^72j. The splitting functions Pap{x/ z) in leading order represent the probability for a parton (3 
with a momentum fraction z of the proton to emit a parton a with momentum fraction x of the parent 
parton. From comparisons of these QCD predictions with data, the quark and gluon momentum 
distributions and the strong coupling constant q;s(M|) can be determined. The equations can be ap- 
plied in 3rd order (NNLO) perturbation theory due to a major progress, achieved recently, in pQCD 
calculations [7^. Most of the parton distribution analyses, as those presented subsequently, are still 
done in NLO. The next order will be important to match the still increasing accuracy of further DIS 
measurements from HI and ZEUS. In the upcoming NNLO analyses, the dependences of the fit results, 
for example, such as that for as(M|) , on the chosen scales are expected to be reduced significantly. 

The HERA data are at sufficiently high at large x such that higher twists, power corrections 
oc Q"^ can be neglected. They extend to such low values of x, however, that issues as to whether 
the evolution equations need to be modified to include Inl/x terms have been under investigation for 
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many years [71]. Besides attempts to replace the DGLAP equations by different equations [5H [75] , 
corresponding to modified rules of multi-gluon emission at low x, the effects have also been investigated 
of their modification with a resummation of small x logarithms [S2]. These studies, so far, have not 
lead to definite conclusions. For the resolution of these questions, an increase in energy beyond HERA, 
as is considered with the LHeC [76] , is very desirable in order to reach even lower x and to perform such 
investigations in a region of large enough where as is small. 

The following sections present a brief overview on the most important pQCD analyses performed so 
far by ZEUS and HI to examine the partonic content of the proton. 

4.2 Determinations of Parton Distributions 

The quark and gluon distributions in the proton define the various structure functions and their x, 
dependence as given in Equation [T71 [31] and [321 The principal method of extracting a complete set 
of quark distributions, q{x,Q^), and the gluon distribution, xg{x,Q^), requires a parameterisation of 
their x dependence at some initial value of = Qq, of a. few GeV^, as input to the DGLAP evolution 
equations and to the factorisation equation (Equation [17]) leading to the structure functions. This 
procedure defines a pQCD prediction of the cross-sections at all and x, which can be confronted 
in a minimisation procedure to the data under study. Two DGLAP evolution programs QCDNUM 
[77] and QCDFIT ^78j, which were independently developed, are used most often in the HERA QCD 
analyses. They have been cross-checked against each other and against a Mellin transformation code 
at both the HERA Physics Workshop in 1995 and subsequently at the HERA-LHC Workshop in 2005. 

The ZEUS and the HI Collaborations so far focused mainly on analyses of their own data, as are 
described subsequently. There have been different attempts for global determinations of the parton 
distributions, such as those by the MRST [79j group, by the CTEQ [80] group and by Alekhine [81], 
in which much wider sets of data, in particular those from lepton-nucleon fixed target experiments and 
the Drell-Yan measurements, are considered. Unlike the fixed-target data, the HERA data constitute a 
complete set of NC and CC cross sections over a very wide range in x and Q^. The additional constraint 
obtained from having both e~^p and e~p data allows the extraction of sets of parton distributions that 
approach the accuracy of global determinations. There are specific regions for which HERA data do 
not have sufficient accuracy; examples are the high x behaviour of the u/d valence quark ratio, which 
is better derived from Ip and ID data, and the u — d difference, which is clearly observed at medium 
X ~ 0.01 [82] in Drell Yan data but was only barely discernable in the initial HI data [83]. This could 
have been accessed at low x with electron-deuteron data at HERA. 

An important consideration in doing pQCD pdf fits is the treatment of systematic errors. In the HI 
fits, the correlated systematic uncertainties are treated as fit parameters; this means the optimum value 
of systematic shifts of the cross-sections are determined through the fit. After the initial fit, the model 
assumptions (the shape of the pdf parameterisation, for example) in the fit are varied to determine the 
"model uncertainties" which are separate from the systematic uncertainties of the data. In the ZEUS 
fits, on the other hand, the initial fit that determines central values of the parameters is made with the 
systematic offsets set to zero. Then an "envelope" of n systematics errors is determined by repeating the 
fits 2n times after off-setting the cross-sections to the +1 standard deviation and -1 standard deviation 
values allowed by the systematic uncertainties and then adding the changes in quadrature B The model 
uncertainties are then found by repeating the central fit with the variation on model uncertainties. The 
uncertainties of the HI pdf results are dominated by the model uncertainties, whereas the systematic 
uncertainties dominate the ZEUS data. It is remarkable that the uncertainties of HI and ZEUS fits 
are, in the end, about the same [M] . 

A further challenge for the QCD pdf determinations consists in the treatment of the heavy flavours. 

•^In practice this is done through a matrix inversion procedure [78j . 
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At HERA, the charm and beauty contributions to the cross section are significant (see below) and the 
data are in regions which include the c and h thresholds. A quark which is heavy, compared to 
the QCD parameter A ~ 0.3 GeV, becomes light at very high compared to M^^. Correspondingly 
there exist various theoretical prescriptions on the treatment of c and 6, as light quarks, as variable in 
their appearance or as heavy quarks, from fusion (see [85] for reviews). Fits have been presented 
with different heavy flavour treatments or versions of the variable flavour scheme (VFNS). ZEUS (Sec- 
tion l4.2.1"]) . has speciflcally used the VFNS scheme of ^86i|- HI has presented a flt to data at lower 
using the flxed flavour scheme, in the determination of as (Section l4.3p . or assumed the heavy quarks to 
be light in the analysis, described in sect. 14. 2.21 which focused on their high data. These and further 
differences lead to somewhat different results. Work is in progress, both by the global flt groups and 
by ZEUS and HI, to further pin down such differences and arrive at flts, desirably at NNLO, based on 
the most plausible assumptions as well as on new and possibly combined data sets. 

The following sections describe briefly the different flts of ZEUS and HI in the determination of the 
parton distributions. In Section the influence of di-jet data is illustrated on the determination of the 
parton distributions, in particular on xg at medium x. Important insight to the behaviour of the gluon 
distribution at low x is expected from the forthcoming direct measurement of the longitudinal structure 
function. 

4.2.1 The ZEUS Approach 

The ZEUS Collaboration performed, in recent years, three major NLO QCD analyses, a ZEUS-S (stan- 
dard) flt [87], which includes the flxed-target DIS data in the flt, a ZEUS-0 (ZEUS data only) flt [88] 
and a flt to the inclusive data together with di-jet data in DIS and photoproduction, ZEUS- Jets flt. 
The ZEUS-0 flt exhibited the power of the HERA data in coming quite close in accuracy and detail 
to the ZEUS-S flt. Here only the ZEUS-0 flt is briefly presented. The effect of adding jet cross sec- 
tion information, (the ZEUS-jets flt) in particular on the gluon determination, is described in the next 
Section. 

The main idea of the ZEUS-0 pdf parameterisation, following a QCD flt tradition, is to decompose 
the quark contributions into sea quarks and valence quarks. Four quark distribution combinations were 
used: 

S* = 2((i -|- U -|- s -|- c), A = — M, dy = d — d, Uv = u — u. (33) 

These distributions (times x) and xg are parameterised as four-parameter polynomials of the type 

xf{x) =pixl{l~xY''{l+p^x). (34) 

The strange sea was coupled to the light sea, x's = O.lxS, and the charm contribution was dynamically 
generated [88]. The most detailed information in the HERA data is on the sea, xS, for which all 4 
parameters were determined, using in addition the momentum sum rule, which flxes the integral over 
xg + xS + xuy + xdv to be one. There is a weak constraint coming from the charged current data, on 
the difference between u and d. Such a difference has been established in the muon-nucleon data as a 
violation of the Gottfried sum-rule and rather accurately measured in Drell-Yan scattering [82] . For x A 
a choice was taken to let only the normalisation parameter pi float and set p2 = 0.5, p^ = PsiS) + 2 and 
Ps = 0. Both valence quark distributions were used flxing p2 = 0.5 and requiring the integral over Uy 
to be 2 and over dy to be 1 (quark counting rule). Finally = was used for xg which constrains the 
high X gluon to be positive. Results from these flts are discussed in detail in [8H|. Subsequently flts were 
performed with additional information on jets. This lead generally to a consistent set of pdf's, both 
with the purely inclusive ZEUS-0 flt, and with the HI flt which is briefly described in the following. 

A comparison of the ZEUS and HI pdf's is presented in Figure [TTl The agreement can be called 
indeed remarkable for it uses different data, different flavour decompositions, different parameterisations, 
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Figure 11: Quark and gluon distributions determined by the ZEUS- Jets and the H1PDF2000 
fits vs. X at = 10 GeV^. Shown are the total error bands which include the experimental 
and the model uncertainties (not separated for ZEUS). One notices the gluon dominance at 
low X at which xg is about 10 times larger than xU, for example. 

different definitions and independent programs imposing differing heavy flavour treatments. More 
detailed inspection reveals that there are differences to be understood, as on the behaviour of xg at low 
X and Q^. 

4.2.2 The HI Approach 

The NLO QCD analysis (H1PDF2000 fit) [H9] of the HI Collaboration of the CC and NC cross sections 
exploits Equations [26] and [28] which suggest a decomposition of the various structure functions into 
the four basic combinations of up {U, U) and down {D,D) quarks and antiquarks (Equation [291) . The 
data used as input are the fixed-target BCDMS data as well as the HI data. The fit was developed 
in conjunction with the publication of the HERAT high data. It was thus decided to use the 
simplest heavy flavour treatment and consider all quarks to be light. As is illustrated in Figure [T^ the 
four effective quark distributions and the gluon distribution xg are rather accurately determined. The 
accuracy achieved so far, at = Ql, for x = 0.01, 0.4 and 0.65, is 1%, 3%, 7% for the sum of up quark 
distributions and 2%, 10%, 30% for the sum of down quark distributions, respectively. It is to be noted 
that in this particular fit the gluon distribution comes out to be rather steeply rising towards lowest x 
at Qq = 4 GeV^, cx a behaviour which in most other analyses is not confirmed. While this low x 

peculiarity had been noted, the fit was still accepted because of a large uncertainty of xg in this region. 
Also, in the higher region for which this analysis was optimised, the fit is in good agreement both 
with previous HI results [67J and with the results obtained in global fits by the MRST |90] and CTEQ 
|91j collaborations (Figure [T2]) . 

The larger x domain is dominated by the valence quarks. In this analysis, these are extracted from 
the differences Uy = U — U and = D — D. Note that this implies the assumption that sea and 
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Figure 12: Determination of the sum of up, anti-up, down and anti-down quark distributions 
and of the gluon distribution in the proton based on the HI neutral and charged current cross 
section data. Left: the parton distributions with their experimental and model uncertainties 
as determined by HI at the starting scale Ql = 4GeV^; Right: for of 10 and 1000 GeV^, 
the HI fit is compared with results from MRST and CTEQ. 



anti-quarks are equal which is contradicted by some non-perturbative QCD models. 

In the HI fit the parton distributions at the initial scale Qq = 4GeV^ are parameterised as xP = 
ApX^^{l — x)*"^ ■ fp{x). The function fp is a polynomial in x which is determined by requiring "x^ 
saturation" of the fits, i.e. starting from fp = 1 additional terms Dpx, Epx'^ etc. are added and only 
considered if they cause a significant improvement in x^, half integer powers were considered in [92]. The 
result for fitting the HI data has been as follows: fg = (l+Dgx), fu = (l+Dux+Fux^), fo = {i+Dnx) 
and fij = f'D = 1- The parton distributions at low x are thus parameterised as xP — > Apx^^ . The 
strange (charm) anti-quark distribution is coupled to the total amount of down (up) anti-quarks as 
s = fcD (c = fcU). In the absence of more detailed information, two assumptions have been made 
on the behaviour of the quark and anti-quark distributions at low x. It has been assumed that quark 
and anti-quark distributions are equal and, moreover, that the sea is flavour symmetric. This implies 
that the slopes B of all four quark distributions are set equal, Bu = B^, = Bjj = Bjj. Moreover, the 
normalisations of up and down quarks are the same, i.e. Ajj{l — fc) = Ajj{l — fg), which ensures that 
d/u 1 as X tends to zero. If the normalisation and slope conditions are relaxed, the uncertainties of 
the pdf's at lower x becomes large [83]. Some constraints can be inferred from fixed target deuteron 
DIS data and will also be derived from asymmetry data at the LHC A high energy collider eD 
experiment would deliver the ultimate answer as to how the sea at low x is composed. 



4.3 Measurements of as in Inclusive DIS 

The scaling violations at large Bjorken x are due to gluon bremsstrahlung. In the valence quark 
approximation the DGLAP equation for F2 takes a non-singlet form, i.e. there is no contribution from 
the gluon to the scaling violations which are solely regulated with as- At low x, on the contrary, the 
dominant process is boson-gluon fusion and the scaling violations are governed by the product of as 
and xg. The HERA data, unlike previous fixed target data at lower energy, therefore have the potential 
of resolving the correlation between as and the gluon distribution for they span such a large range in 
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X and Q^. This has been exploited by ZEUS and HI to determine both the strong couphng constant 
and xg from the inclusive DIS data using NLO perturbative QCD. 

The ZEUS Collaboration has used its standard (S) fit and parameterisations as described above. In 
addition to the ZEUS data, a number of fixed target experiment data was used in addition (see [HH]), in 
order to have better accuracy for the high x region and constraints on the up-down quark asymmetry 
at medium x. The heavy quarks were treated in a variable flavour scheme. Standard cuts were made 
on the input data with > 3 GeV^ and > 20 GeV^ in order to be in the DIS region and avoid 
higher twist contributions, from low and large x, respectively. The input scale at which the pdfs are 
parameterised was chosen to be Ql = 7GeV^. The result obtained is 

a,(M|) = 0.1166 ± 0.0049 (exp) ± 0.0018 (model) (35) 

where the experimental error comprises three uncertainties, 0.0008 from uncorrelated error sources, 
0.0032 from correlated error sources and 0.0036 from normalisation uncertainties. The model uncer- 
tainties stem from variations of the fit parameters, Qmin^ ^min and VF^m; variations of quark 
and gluon distribution function parameterisations. The result changes by 0.001 if the heavy flavours 
are treated in the fixed flavour scheme. Special care was devoted to the effect of higher twists as the 
fit included fixed target data at rather small but large x. With an x dependent parameterisation of 
power corrections oc Q~^, it was observed that as{M^) may be lowered by 0.0032. 

The HI Collaboration developed a dedicated procedure [67] for the determination of which 
differed in several respects from their pdf fit described above. The input data were restricted to HI and 
BCDMS lepton-proton scattering data when it was found that the NMC data, extending to low of 
a few GeV^ introduced a significant dependence of the as{Mz) result on Q^j„ (see also [93]). Since this 
requires the movement of the minimum for the fixed target data up to about 10 GeV^, the accuracy 
of the BCDMS data is so overwhelming that no further fixed target data has to be included, apart from 
consistency checks [67]. The analysis also excluded the low y data of BCDMS from the analysis for 
systematic reasons. The second feature of the HI fit was the attempt to reduce the number of quark 
distributions to be fitted from the usual three to two. This was possible as the fit was essentially on F2 
only, which can be written as a sum over singlet and non-singlet contributions as 

F2 = --x^ + --x/\ = --xV + — -xA. (36) 

9 3 3 9 ^ ' 

The V, A functions are given as 

V = ^[3(f/ + U)- 2{D + D)] = ^u, - + 3(d + s) (37) 

A=^[2{D + D)-{U + U)] = -^u, + ^4 -^u + (d + s) (38) 

for 3 fiavours as heavy fiavours were generated according to boson-gluon fusion. In the simplified 
approximation, Uy = 2dy and u = d = 2s, one finds that V and A are valence and sea quark dominated 
functions, respectively, as in this approximation V = 3uy/2 and A = u. Despite the reduction in the 
pdf space, the fit could impose the usual momentum and quark counting rules and also a non-zero u — d 
difference. The result obtained is 

a,(M|)= 0.1150 ± 0.0017 (exj9) + (morfe/). (39) 

The first error represents the experimental uncertainty of the data sets. The second error includes all 
uncertainties associated with the construction of the QCD model for the measured cross section. A 
number to be noted may be that a 100 MeV uncertainty on the value of the charm mass results in an 
a,(M|) uncertainty of 0.0005. 
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The experimental results of ZEUS and HI are consistent, and with more accurate and the combined 
data sets a further improved measurement will be reached. A rather large theoretical uncertainty of 
the NLO analysis, however, results from the choices of the renormalisation scale /i^ = ■ Q^, and of 
the factorisation scale fij = ruf ■ Q'^ which leads to scale dependent parton distributions. In the MS 
renormalisation scheme both scales are set equal to Q^, i.e. = rrif = 1. The effect of both scales on 
a^was estimated by both ZEUS and HI through variations of the scale factors and ruf. 

The ZEUS Collaboration [94] has varied rrirj by factors between 0.5 and 2. This caused shifts in 
a;s(M|)of ~ ±0.004. Variations between 0.25 and 4, as are often common, resulted in unacceptable 

variation effects and were thus not considered. Similar results were obtained by HI [67] as are 
summarized in Table [3l 

Table 3: Dependence of as(M|) in the HI fit on the renormalisation and factorisation scales 
ruf and m^, respectively, expressed as the difference of as(M|) obtained for scales different 
from one and the central value of ^^(Mf ) =0.1150. The combination m/ = 4 and rrir = 0.25 
is abandoned since the splitting function term cx In [mr/mf)'^ becomes negative at low 
which causes a huge increase of x^- 





rrir = 0.25 


rrir = 1 


rrir = 4 


ruf = 0.25 
mj = 1 
nif = 4 


-0.0038 
-0.0055 


-0.0001 
+0.0005 


+0.0043 
+0.0047 
+0.0063 



The numerical size of the effects depends on the data and Q^in^ the HI and ZEUS observations are 
thus found to be consistent. In agreement with previous studies [93j it was observed that the renor- 
malisation scale causes a much larger uncertainty on a;s(M|) than the factorisation scale. Depending 
on which set of rrir and is chosen, the obtained typically differs by a number much larger than 
1. This suggests that the assumed variations of the scales are too large. It has to be stressed again 
that while there is an arbitrariness in the scale choices, which is inherent to the renormalisation group 
equation, the actual prescription of variations by typically 1/4 to 4 times the scale is ad hoc and needs 
to be revisited. The scale dependence gets diminished in NNLO [73], yet future determinations of 
q;s(M|) from HERA data and possibly from new colliders [95] will push the experimental accuracy to a 
per cent or per mille, respectively. In order to take advantage of such small experimental uncertainties, 
a new understanding of theoretical uncertainties-particularly that of the scale uncertainties-will be 
necessary. 
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5 Jet Measurements 

5.1 Theoretical Considerations 

According to the QCD improved parton model, a DIS differential cross section, da can be written as, 

daep = E / dxfa{x, as) ■ daa{xP, a^, hr, hf)- (40) 

This has the form of a convolution of a partonic, or hard, cross-section da with the proton PDF fa- It 
should be noted that both da and fa depend on as which depends on ^r. 

In order to produce theoretical predictions in practice, semi-analytical programmes such as DIS- 
ENT [96j MEPJET |j97j or DISASTER [98] are used. The partonic cross-section, da, above a specified 
scale hr, is calculated to NLO in as by these programs. The PDF of the proton, fa, is normally taken 
from the results of global fits such as the MRST [HS] or CTEQ [lUU] sets of PDFs. These global sets 
do not include the HERA jet data in their fits. These PDF fits are made using a fixed value of as, but 
PDFs are, in principle, a function of as, as indicated in Equation l40l 

In order to produce jet cross-sections from the partonic cross- sect ions produced by DISENT and 
DISASTER, a jet algorithm0, the inclusive kx algorithm |1031I104] in case of HERA jet analyses, are run 
on the final partons to produce the "parton-level" jet predictions to NLO. The actual measurements, 
however, are not made at the parton-level, so additional corrections are necessary. 

Monte Carlo(MC) programs currently available for DIS event generation, such as LEPTO [105] or 
ARIADNE [106(1107] . only contain the LO matrix elements. They are, however, supplemented by ad-hoc 
additions of NLO processes as well as a simulation of higher order processes through the Parton-Shower 
model. In addition, the partons are "hadronised" using the LUND [108] model. Thus, after running 
the jet algorithm on the events produced by MC programs, both "parton-level" and "hadron-level" MC 
predictions can be obtained. 

In order to obtain the "hadron-level" prediction to order a^, a multiplicative "hadronisation cor- 
rection" is applied to the parton-level prediction. In most cases the hadronisation correction is the 
ratio of hadron-level to parton-level MC predictions. While there is no rigorous justification for this 
procedure, it is a reasonable one if a) the parton distributions in the MC are close to those of the fixed 
order calculation, b) the MC reproduces the experimental measurements reasonably well and c) the 
hadronisation correction is small (5-10%). Normally at least two MC programs are compared and used 
to evaluate a systematic uncertainty due to this procedure. 

It is convenient to symbolically write the above procedure to obtain the jet cross sections in DIS as: 

^^ef = E / dxfa{x, as) ■ dal''\xP, asifJ^n), fiR, fip) ■ (1 + Shad), (41) 

where rfa^"^* is now the "jet" , rather than the partonic, hard cross-section and 6had is the hadronisation 
correction. 

In the case of photoproduction, where ~ (i.e. an almost real photon collides with the proton), 
the photon PDF, as well as the proton PDF is factorized. Thus, the photoproduction (7p) jet cross- 
section can be written as: 

do'i/p = EE/ j dxpfa{xp,^l],\as)■dx^fb{x^,^],]as)■dal''\xp,x^,as{^lR),^^R,^^F)■{^ + 5had), 

a=q,q,g b=q,q,g 

(42) 

*Jet algorithms produce a small number of "jets" from a larger number of final state objects. The four-momenta of 
jets should correspond closely to those of hard partons. Discussions of theory and application of jet algorithms at collider 
experiments may be found in [1011 1102j . 
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Figure 13: Measurement of jet cross-sections at HERA. The filled dots are the data. The 
inner error bars represent the statistical uncertainty. The outer error bars show the statis- 
tical and systematic uncertainties in quadrature. The NLO QCD calculation, corrected for 
hadronisation effects (and exchange when the calculations do not take electroweak effects 
into account) is shown. Left-top) The differential cross-section da/dQ"^ for the inclusive jet 
production. Left-bottom) the ratio between the measured da jdQ^ and the NLO QCD cal- 
culation. The shaded band displays the total theoretical uncertainty. Right) The differential 
cross-section da /dE^^^^^ in bins of Q^, where E^^^et is the transverse energy of the measured 
jet in the Breit frame |117j . 



where fb is the PDF of the photon, is the Bjorken x of a parton of the photon, and the convolution 
extends over partons a in the proton and h in the photon. 

At LO (i.e. to zero-th order in a^), the photoproduction process can be thought of as being in 
two pieces; the "direct" process in which the photon couples directly to the proton (x^ = 1), and the 
"resolved processes" where a parton from the photon interacts with the proton (x^ < 1). While this 
distinction does not have a strict meaning at higher orders, it is convenient to distinguish these two 
processes defined by means of a cut on x^, usually at 0.75 or 0.8 by convention. The "direct" process 
loosely corresponds to the case where = 1, and is not very sensitive to the hadronic nature of the 
photon through its PDFs. 

The NLO theoretical prediction for photoproduction is obtained in much the same way as in the 
case for DIS. The NLO analytical programme often used is that of Frixione and Ridolfi |109] . The most 
usual LO MC programs are PYTHIA [LlOl and HERWIG [TTT] . 

Jet production in the region of transition between photoproduction and DIS has also been extensively 
studied at HERA. The reader is referred to the following HERA results for further information on this 
subject \112\ . 
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Figure 14: The measured photoproduction dijet cross-sections differential in ETjeti, the 
transverse energy of the highest energy jet, is shown separately for > 0.75 (left) and 
< 0.75 (right). The NLO prediction corrected for hadronisation effects is shown calculated 
using the GRV-HO and CTEQ5M1 PDFs for the photon and the proton, respectively. The 
hatched band represents the theoretical uncertainties. The prediction using an alternative 
photon PDF, AFG-HO [m], is also shown. 



5.2 Jet Cross-Section Measurements 

The experimental jet measurements are made at the "detector-level". This means that the energy 
deposits in the calorimeters and tracks found in tracking detectors are subjected to the jet algorithms 
to determine the jets and their four momenta. The jet algorithm used is the same as that used in the 
theoretical predictions described in the previous section, i.e. the inclusive kjn algorithm. 

A Monte Carlo simulation is used to correct the measured detector-level jet cross-sections to the 
hadron-level cross-sections, which can be compared to the theoretical predictions. Events generated by 
MC generators such as LEPTO or ARIADNE for DIS, and PYTHIA and HERWIG for photoproduction, 
are processed through detector simulations based on GEANT [114J to give a sample of fully simulated 
events. The correction factors are calculated using the difference between the detector- and hadron-level 
event rates of the fully simulated event samples. The event selection and kinematic cuts are devised in 
such a way to keep these corrections relatively small. Typically the corrections are < 30%. 

Figure [13] shows the measurement of inclusive jet cross-section in DIS [1151 1116j . in this 
function of Q^. The uncertainties of the measurement, except at the highest Q^, where statistical errors 
are large, are dominated by the jet energy scale uncertainties. The energy uncertainties are 1 to 2% at 
HERA experiments which translates to 5 to 10% uncertainties in the cross-sections. 

The measured cross-section in Figure [13]^ right), for example, is compared to the predictions of NLO 
QCD, in this case based on the DISENT programme. The inputs to the calculation are as for the 
partonic cross-section, chosen to be 0.1175, the factorisation and renormalisation scales (set to be either 
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Figure 15: Left: The inclusive trijet cross-section measured as a function of Q^. The pre- 
diction of pQCD in leading order (dotted line) and in NLO (solid line) and without hadro- 
nisation corrections (dashed line) are compared to the data. Also shown is the ratio of the 
measurement and the theoretical prediction. Right: a) The inclusive trijet cross-sections 
as a function of Ei^^^ with the jets ordered in Ej^^^. The cross sections of the second and 
third jet were scaled by factors hown inthe figure for readability. The NLO prediction uses 
CTEQ6 PDFs of the protons. The other details are the same as in Figure [131 b) c) and d) 
are ratios of the measurements to the NLO predictions. 



Q or Ej^ ) and the proton PDF, chosen to be the version of MRST99 |118] extracted with set to 
0.1175. The bottom plot shows the hadronisation correction, or parton-to-hadron corrections, which 
converts the parton-level predictions of DISENT to the hadron-level predictions to be compared to 
data. This particular correction was determined using ARIADNE and LEPTO MC programmes. The 
comparisons to the theory are generally satisfactory, and the residual difference can be attributed to 
the input parameters of theory, namely and the PDF of the proton. 

Very similar considerations apply to the study of photoproduced jets at HERA. One difference is 
that the PDF of the photon is an additional input to the theoretical prediction as discussed above. 
Figure m] shows the comparison of the theory and data |119] for photoproduced dijets for the "direct" 
and "resolved" samples. Many calculations of jet photoproduction at NLO exist, all of which have been 
compared to each other and agree to within 5-10%. The calculations shown here are due to Frixione 
and Ridolfi [109]. The photon PDF used in this case is GRV-HO [T20lfT2T] . 

The DIS inclusive and photoproduction dijet measurements discussed in some detail above are two 
of the most basic jet measurements at HERA. The follow-up measurements to these can be found in 
[122^ 1116^ 1123^ 1124] . An example of a more complex measurement is that of multijet in DIS which is 
shown in Figure [15] [125^ 1126] . The NLO calculations describe the data satisfactorily. There exist many 
more detailed measurements which are listed here [1271 11281 1129] for reference. 

In addition to the measurements of cross-sections, there are tests of pQCD that can be performed 
at HERA using jet measurements. An example is the angular distribution of photoproduced dijets. 
Figure [16] shows a measurement by the HI Collaboration ^122j of the photoproduction dijet cross- 
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Figure 16: Photoproduction dijet cross-section as a function of |cos^*| for events with the 
dijet mass greater than 65 GeV. The QCD NLO calculations with (solid line) and without 
(dashed line) hadronisation corrections, 6had and for the PYTHIA Monte Carlo (dotted ine) 
scaled by a factor 1.2 are shown. The inner (hatched) band is the scale uncertainty the outer 
(shaded) band is the total uncertainty. The cross-section are shown for the two regions of 
enhancing the resolved (left) and direct (right) photon contribution. 



section differential in cosine of the angle of dijets, | cos6'*|, in their centre-of-mass, with respect to the 
beams. The difference of steepness of the | cos 9* \ dependence between the cross- sect ions from resolved 
and direct samples is predicted by the dominance of gluon exchange in the resolved process and is 
related to the spin of gluons and quarks. 

5.3 Tests of pQCD and Determination of as 

The strong coupling constant is extracted by a fit of the theoretical predictions to the data with 
as as a parameter of the fit. At HERA, the dependence of PDFs on may be handled consistently 
by a use of sets of PDFs determined at different as such as the MRST99 set. The best fit is found 
by simultaneously varying as in the hard cross-section as well as the PDFs. The PDFs determined at 
discrete values of as are used to interpolate to the intermediate values. In this way, as is obtained in a 
consistent manner taking the correlations into account. 

The extraction of as can be done in bins of Et (or Q'^) to show its running as a function of the 
renormalisation scale. This is shown in Figure fTW left) for a variety of jet measurements at HERA |129j . 
In all measurements of as from jets at HERA, theoretical uncertainties dominate the systematic error, 
the largest of these are due to unknown higher order contributions, evaluated by a variation of the 
renormalisation scale, fipt- 

A compilation of the ag values has been made; a combined HERA as has been derived |13Uj and is 
shown in Figure [T71 (right). The overall experimental uncertainty of the combined measurement is 1%, 
the theoretical uncertainty, which is dominated by the scale uncertainty of the NLO theory, is about 
5%. As was noted in the inclusive case, there is no firm procedure which would define the theoretical 
uncertainty which remains ad- hoc to a large extent. This combined HERA measurement has been 
used in a determination of the world average of the strong coupling constant }132] of 0.1189 ± 0.0010. 



Recently, a determination of the strong coupling constant was made using a simultaneous fit of the 
HERA data from [HS] and to give a value of 0.1198 ± 0.0019(exp) ± 0.0026(th) [ISS]. In 

this analysis only data at large scales were used. This halved the theoretical error but enlarged the 
experimental one. Whether theoretical uncertainties can be settled in a more satisfactory way remains 
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Figure 17: Left: A compilation of asiE^jet, Q^) measurements at HERA determined from the 
QCD fit of tlie measured da/d^Exjet, Q"^) as a function of Et^jet or Q"^. The band shows the 
current world average (excluding HERA data). Right: The combined HERA measurements 
from |13Uj . The dashed line and the band shows the result of a global analysis (not using 
the HERA data) from [m]. 



to be seen. 



5.4 Simultaneous Determination of PDF and ag using Jet Data 

An examination of Equations HQ] and Hi] shows that the jet cross-sections at HERA are sensitive to 
both proton PDFs and Ug- In the previous section, the dependence on proton PDFs was treated as a 
systematic uncertainty in determining from the jet data. One may also think of determining the 
PDFs from jet data using determined elsewhere. Indeed, the simultaneous determination of proton 
PDF and in a combined NLO QCD fit of DIS inclusive and jet cross-sections (and possibly other 
data) is the most consistent approach. 

The inclusion of HERA jet data, or indeed any jet data in an NLO QCD fit, such as the one described 
in Section 4, is not conceptually difficult; the problem is technical. While the DGLAP equation can be 
rapidly evaluated in a programme such as QCDNUM [77j, the NLO jet cross-section calculations are 
time consuming enough to make using them in an iterative fitting procedure impractical. Only recently 
have technical difficulties been overcome, making this procedure possible. 

In 1 134] , the hard jet cross-section, da^^^ in Equation |lT]is pre-calculated on a grid of four dimensions, 
X, the proton momentum fraction, a^, Aii?, the renormalisation scale and ^p-, the factorisation scale. By 
an appropriate choice of the grid points at which the pre-calculations are carried out, an accuracy of 
better than 0.5% with respect to a full calculation can be achieved using interpolations for all relevant 
cross-sections. In this way, the full NLO jet calculations can be incorporated into a NLO QCD fit on 
the same rigorous footing as the DGLAP equations for inclusive cross-sections. 

The results of the fit using the ZEUS DIS inclusive cross-sections as well as both DIS inclusive 
jet cross-sections and the direct photoproduction dijet cross-section are shown in Figure [18] In this 
instance, the fit is made with the value of fixed to the world average value. There is a significant 
improvement to the PDF uncertainties, particularly at values of x above 10~^, which is expected since 
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jet production is dominated by gluon-photon fusion process in this range of x at HERA. 
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Figure 18: Left: The gluon and sea PDFs extracted from the ZEUS- JETS fit. The inner 
cross-hatched error band shows the statistical and uncorrelated systematic uncertainty, the 
grey error bands show the total uncertainty including experimental uncorrelated and cor- 
related systematic uncertainties, normalisation and model uncertainty. Right: The total 
experimental uncertainty on the gluon PDF for the ZEUS- JETS fit (central error bands) 
compared to the total experimental uncertainty on the gluon PDF for a fit not including the 
jet data (outer error bands). The uncertainties are shown as fractional differences from the 
central values of the fits for various values of Q^. 



The next step is to simultaneously fit as as well as the PDFs. In a QCD fit using only the structure 
function data from HERA, is so far weakly constrained as is shown in Figure dni The same figure 
shows the strong constraint on jet data add to the fit. The uncertainties of the PDFs also remain 
stable and relatively small for the simultaneous fit as shown in the same figure. The extracted value 
of as{Mz) = 0.1183 ± 0.0028(exp) ± 0.0008(model), with an additional uncertainty of ±0.005 due to 
unknown higher order effects to the NLO theoretical calculations, is very close to the value used in 
the fixed-a^ fit, and thus does not change the central values of the PDFs significantly. The level of 
the precision of this determination of is as good as any existing single measurement; it is also, 
so far, the only accurate one that is derived from HERA data only, with minimal assumptions in 
the PDF fit parametrisation that derive from outside data. The other jet measurements use PDF fit 
parametrisations derived using world data. 
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Figure 19: Left: The profile as a function of a^iMz) for the fit to ZEUS data with (black 
dots) and without (clear dots) including the jet data. The ordinate is given in terms of the 
difference between the total and the minimum for each fit. Right: The gluon distri- 
butions extracted from the simultaneous fit to ZEUS DIS and jet data. The uncertainties 
on these distribution are shown beneath each distribution as fractional differences from the 
central value. The outer error band shows the total uncertainty and the inner band the 
uncertainty in case is fixed in the fit. 



5.5 Summary and Outlook of Jet Measurements at HERA 

The jet cross-section measurements at HERA have achieved about 5% precision, with the main contri- 
bution coming from the understanding of the jet-energy scale at the level of 1-2%. With this precision, 
a meaningful comparison can be made to the QCD NLO predictions of jet predictions. Indeed, in most 
cases, the precision of the experimental measurements are such that the theoretical part dominates the 
uncertainties. 

Based on the good description of the data, QCD parameters can be extracted. The strong coupling 
constant is determined to about 1% experimental uncertainty (5% theoretical uncertainty) from the 
current measurements. The jet measurements help to constrain the gluon PDF when combined in a 
simultaneous fit with the DIS inclusive cross-sections. The current measurements constrain the gluon 
to about ±10% at x around 0.05 and around 200 GeV^, based so far on the ZEUS jet data only. 

All of the results discussed here come from the HERA I data taking phase which comprises about 
120 pb^^ per experiment. While many of the jet measurements are now limited by the theoretical 
uncertainties, there are areas, such as the measurements at the highest which will benefit from the 
800 pb~ of which is available from the HERA II data taking period. The projections are that this data 
should be able to improve the uncertainty of the gluon PDF in the x range of 0.1 to 0.05 by about a 
factor of two [135j . 
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6 Heavy Flavours 

6.1 Introduction 

Processes involving heavy quarks, with mass, Mq ^ ^qcd, are, in principle, amenable to pQCD 
calculations. At HERA, the production of charm and beauty particles are measured and compared 
to theoretical predictions. While Mq is one QCD scale, there are two other relevant scales in heavy 
quark production processes at HERA. These are the virtuality, Q^, of the exchanged photon, and 
the transverse energy, Et, of the final state particles or jets. The relative size of these three scales 
is important in choosing what kind of calculation is appropriate for comparison to which data. For 
example, in the limit, ^ Mq, an approximation where Mq = may be appropriate. On the other 
hand, since there are no quantitative predictions of what "much greater than" may mean, the study 
of heavy quark production, as is often the case with QCD, is an iterative one where the appropriate 
theory is arrived at gradually as better and better data are compared to the theory. 

6.2 Theoretical Calculations of Heavy Quark Production 

In the conventional co-linear approximation, heavy quark production in DIS cross-section can be written 
generically as 



If the renormalisation scale, hr, usually taken to be the virtuality Q^, is ^ Mq, then the heavy quark 
should be treated as a massless quark on the same footing as the other quarks in the PDF in the 
framework of In/x^ resummation (DGLAP approximation). In this case, fa is the PDF of the heavy 
quark, and the cross section da is the photon- (photoproduction) or electron- (DIS) quark cross-section 
for a massless quark. This type of calculation is called "massless" . 

If the scale, ~ Mq, then it is appropriate to treat the heavy quark separately as a massive quark. 
In this case, the sum over a in the above equation will be for all partons except the heavy quark. The 
cross-section da, will be for a production of a massive quark from an electron colliding with a (massless) 
parton at a fixed order in a^. This type of calculation is called "massive". 

There are three types of treatments of heavy quark production in DIS based on co-linear factorisation. 

• Zero-Mass Variable Flavour Number Scheme (ZM-VFNS), is a purely "massless" treatment where 
the heavy quarks are treated as massless quarks which become "active" in the proton PDF at 
the production thresholds of a QQ pair. This is the oldest and still the most common theoretical 
treatment. As discussed above, this is, in principle, incorrect near the production thresholds. 

• Fixed Flavour Number Scheme (FFNS) is a purely massive treatment. This is, in principle, 
incorrect at high Q^. The most common treatment uses the analytic programme HVQDIS [136j . 

• Matching Schemes. Sometimes these are called Massive Variable Flavour Number Scheme (Massive- 
VFNS) or Generalised Mass Variable Flavour Number Scheme (GMVFNS). There are several pre- 
dictions which match the massless and massive schemes at an intermediate value of Q^. Examples 



of these are Roberts-Thorne (RT) [1371 [I38], ACOT [I39], Kniehl et al. [HO] and ResBos-HQ [Til] . 



The cross-sections, calculated to order a^, are available in all three schemes. Some of the predictions, 
such as those of HVQDIS and ResBos-HQ, are of semi-analytical type similar to those discussed in 
Chapter 4 for jet production. In these detailed prediction for comparison to the data of 

measurements of production of hadrons containing heavy quarks is possible. Such predictions are made 
using a fragmentation function |142l I143[ 1144] that give the probability of producing a particular hadron 
given an initial heavy quark. Also, a cross-section of jets associated with heavy quark production can 
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Figure 20: Differential D* ep cross-sections as a function of a) Q^, b) x, c) pt{D*) and 
d) ri{D*) compared to calculation of HVQDIS CASCADE and a Massive- VFNS calculation 
from Kniehl et al. [150j. The inner error bars indicate the statistical uncertainty, while the 
outer show the statistical and systematic uncertainties in quadrature. Figures a), b) and c) 
also present the ratios R = o'theory/o'data, taking into account their theoretical uncertainties. 
See [HH] for details. 



be predicted in the same way as for ordinary jet production discussed in Section 5. Many predictions 
give only the x and dependence of the charm cross-section, and thus can only be compared after 
some, more or less large, extrapolation of the data (see also the end of this Section). 

The theoretical treatment of the photoproduction of heavy flavours is similar to that of DIS but 
complicated by the second "hadron" , the real photon. In this case the photon-proton cross-section can 
be represented as; 

da = J2j J dxpdx^faixp, ■ fb{x^, /x^; a,) ■ d&HQixP, /in, I^f, Mq), (44) 

a,b 

where fb is the PDF of the photon. In practice, the measurements use jets associated with heavy quark 
production to determine the kinematics-in particular the Bjorken variable of the photon, x^. 

As in the case of DIS, predictions in FFNS, ZM-VFNS and Massive- VFNS exist for the photopro- 
duction case. The renormalisation scale /ir is usually chosen to be some combination of the quark mass 
Mq and the transverse momentum of the produced heavy quarks, . 
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Figure 21: The measured F!^'^ in bins of as a function of x. The data are shown with 
statistical uncertainties (inner bars) and statistical and systematic uncertainties added in 
quadrature (outer bars). The lower and upper curves show the fit uncertainty propagated 
from the experimental uncertainties of the fitted data. 



The predictions discussed above are all based on co-linear factorisation. A different approach is to 
use kt factorisation and the CCFM equations ^145j , which, unlike the DGLAP equations, resum terms 
proportional to In 1/x as well as to In The actual predictions in this approach are made by the Monte 
Carlo program CASCADE |146j that uses the off-shell photon-gluon fusion process convoluted with the 
gluon distribution unintegrated in kt, obtained via a fit to HERA F2 data [147j . These predictions are 
currently available only at order as- Like the predictions of HVQDIS, the CASCADE programme is 
capable of providing predictions at the parton level, and is therefore easily comparable to measurements. 

The reader is referred to jl48J for a more complete discussion of the theory of heavy quark production 
than is possible here. 

6.3 L>* Cross Sections in DIS 

While many different kinds of charmed hadrons have been measured at HERA, the measurements of 
the D* mesons have, by far, the highest statistical significance. This is due to the well-known ease of 
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Figure 22: Inclusive D* cross-sections as a function of pt{D*) and ri{D*) compared to the 
QCD NLO calculations of FMNR and GMVFNS. For FMNR the beauty contribution is 
shown separately for pt{D*). The ratio of the measurement and the prediction, R, is also 
shown. 



reconstructing the D* meson, using tracking detectors alone, by exploiting the "slow" pion from the 
decay, D* ^ D + 11^. 

Figure [20] shows a recent measurement of the D* cross-sections in DIS [149]. The theoretical predic- 
tions, which all describe the data reasonably well are of the FFNS type, HVQDIS, as well as the CCFM 
based CASCADE and a massive VFNS calculation. The predictions use a range of proton structure 
functions and charm quark masses. It is interesting to note that the description of FFNS appears to be 
somewhat better than the others even at rather high values near a hundred GeV^ despite the fact 
that FFNS should be, in principle, incorrect when ^ Mq. Indeed this good description is known to 
hold up to several hundred GeV^ [151] . 

6.4 Derivation of the Charm Structure Function 

The measured charm hadron cross-section can be extrapolated to the full charm cross-section differential 
only in x and Q^. This, in turn, can be interpreted as the charm part of the cross-section, described 
essentially by the structure function F2^. Thus, 

^ = + (1 - yfWfix, - y-Ft{x^ Q^)}. (45) 

Since the D* cross-section is measured in a kinematic range limited in and ?], a method of extrapola- 
tion is needed. A most straightforward way to determine this is, for the cross-section cTj(ep D*X) 
in the i-th bin of x and Q^: 



pec ( q2^ _ (^i,meas{ep ^ D*X) ( ^ . O'^) (46) 
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Figure 23: Charm dijet cross-section da/dx° compared to QCD predictions in the massive 
scheme, with (sohd hne) and without (dotted hne) hadronisation corrections apphed. The 
theoretical uncertainties (hatched band) come from the change in scales simultaneously with 
the change in charm mass. 



where Xi and Qj are chosen to be within the measured bin. The uncertainties in this extrapolation are 
determined by varying the parameters of the theoretical models being used. While these uncertainties, at 
10-20%, are usually smaller than the uncertainties already present in the D* cross-section measurement, 
it is worth keeping in mind that the extrapolation factors range from about 1.5 up to nearly 5. Thus 
reliance on the correctness of the theoretical prediction of the cross-section in rj and Pt is relatively 
large. This is overcome with impact parameter based measurements as are discussed below. 

Figure [21] shows a compilation of the measurements of F2'^ based on measurements of D* mesons |15H 
I152j . The measurements are compared to a prediction of a PDF obtained from a fit using a VFNS 
scheme. It should be noted that the comparison of theory to data in the form of F^*^ derived from the 
D* measurements do not contain more information than the comparison of the differential cross-sections 
in the measured kinematic range to the theory of the kind shown in Figure [201 

6.5 D* Production in Photoproduction 

A recent measurement |153j of charmed meson production in photoproduction is shown in Figure [221 
The theoretical predictions of the massive scheme (FMNR |154j ) have a tendency to have too shallow 
a slope in pj- of the D* and be below the data at forward rj of the D*. The predictions of the variable 
flavour number scheme (GMVFNS) [155] . while describing the data, have a very large uncertainty 
associated with effects of missing higher orders. 

In photoproduction, one factorizes the photon structure as well as the proton structure in order to 
apply pQCD. Therefore, in order to compare the predictions to the data in more detail than above, it 
is necessary to measure, in addition to the total 7p cms energy, W, and the kinematics of the charmed 
hadron, quantities that can determine the kinematics of the partons in the photon. This can be achieved 
by measuring the jets produced in the events with charm hadrons. The jet finding is done in the same 
way as is described in the previous section. The jets found may contain the charmed hadron produced 
in the event. 

Figure [231 shows a measurement of dijet photoproduction associated with a D* [156]. The cross- 
section is plotted differentially in x^, the Bjorken variable of the photon. The theoretical prediction, a 
massive calculation, tends to fall below the data in the region of low x^ indicating that the description 
in the resolved region is not adequate. Figure [2^ shows the same data now plotted separately for high x^ 
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Figure 24: Cross-section da/dA(f)^^ separated into resolved enriched {x°^'^ < 0.75) and direct 
enriched {x°^'^ > 0.75) region. The theoretical uncertainties come from the change in scales 
simultaneously with the change in charm mass. The beauty component is also shown (lower 
histogram) . 



(direct component) and low (resolved component) in terms of A4>jj, the azimuthal separation of the 
two jets. It is seen that the theory tends to fall below the data away from Acpjj = tt (i.e. back-to-back 
jets) particularly in the resolved region. Since A0jj is equal to tt to LO in as, NLO is the order at which 
the first non-trivial prediction is given. A failure to describe this distribution may be an indication of 
a need for higher order terms in the calculation. 

The comparison of the data on photoproduction of charm with the present theoretical predictions 
indicate some inadequacies of the theory. The indications, both from relatively large renormalisation 
uncertainties and from examination of differential distributions, are that the theory still requires the 
inclusion of higher order terms in order to give a good description of the available data. 



6.6 Photoproduction of Beauty Particles 

The measurements of beauty photoproduction are based on the semi-leptonic decay of the B-hadron. 
Because of the large mass of the 6-quark, the leptons from the decay tend to be produced with large 
transverse momentum with respect to the momentum of the rest of the associated jet. In this case, the 
signal is extracted from the distribution of momenta of the leptons relative to that of the jets. Alterna- 
tively, a D* may be explicitly reconstructed with a decay lepton and their momentum separations. In 
the latter case, which exploits the specific decay chain B — >■ D*^i>, further constraints, such as charge 
of the D* and the fi, and the mass of the pair, can be used in the selection. 

Figure [25] is a compilation of all currently available measurements |157i I158i I159[ I160[ I16H 1162] 
of 6-photoproduction, measured via its decay, converted to a 6-quark cross-section differential in the 
transverse momentum of the 6-quark. It can be seen that while the massive calculation gives a predic- 
tion which is somewhat lower than the measurements, particularly at lower ranges of p^, the general 
description is good. It can also be seen that the kt factorisation calculations from |163j is very similar 
to that of FMNR. 
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Figure 25: A summary of all available HERA data on beauty pliotoproduction plotted 
as the cross-section da/dp\^ and compared to NLO QCD in the massive scheme and a Ut 
factorisation calculation. 



6.7 Measurements of Heavy Quark Production using Vertex Separation 

Qualitatively different measurements of heavy quark production at HERA than discussed above are 
achieved by the use of the precision tracking made possible by the use of silicon vertex detectors. These 
devices have been installed as upgrades to the HERA detectors. The HI vertex detector |164] was 
installed in 1996-7 and the ZEUS one \^Ml in 2001. 

The spatial impact parameter resolution of the silicon detectors, in the order of tens of microns, 
makes possible the use of the long lifetimes of c and h flavoured hadrons for selection of these events 
without explicit reconstruction of the hadrons. The events containing heavy quarks are distinguished 
from those containing only light quarks by reconstructing the displacement of tracks from the primary 
vertex. The major advantages of this method compared to that of reconstructing particular c or 6 
hadrons are that a) there are little or no uncertainties associated with properties of particular hadrons, 
such as its fragmentation function or decay branching ratios, and b) the severe limitation in the mea- 
surable kinematic range in and r] due to the use of central tracking chambers for the reconstruction 
of hadrons is removed, the acceptance for the silicon detectors is close to 100%. The latter point is 
very important when extracting the total c and h cross-sections (and therefore F2'^ and F2^) since the 
extrapolation factors from previous analyses can be as large as 4 or 5 and, thus rely heavily on the 
correctness of the models. 

The first such type of measurements of the heavy quark structure functions |166l I167j was performed 
by HI and is shown in Figure [261 together with previous HI and ZEUS measurements of F^^^ that rely 
on the reconstruction of D* mesons. It is notable that the F^*^ measurements from D* reconstruction 
agree well with the newer results, which is an indication of the reliability of the extrapolations used in 
the older results. While there is an earlier measurement of a{hb) in DIS based on the relative transverse 
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momenta of muons |168] this is the first measurement at HERA of in a relatively wide kinematic 
range. 
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Figure 26: The measured heavy quark structure functions Fg'^ and shown as a function of 
for various x bins. The solid points are the measurements utilising the silicon vertex de- 
tector of HI. The open points on F^*^ are older measurements made using the reconstruction 
oi D* mesons. 

Both ^2^^ and rise steeply as a function of at fixed x, especially at small x, reflecting the 
large gluon density there. The predictions of NLO QCD in the Massive VFN schemes describe the data 
well, except at the lowest x for Fg'^^ where the predictions tend to undershoot the data. 

6.8 Summary of Heavy Quark Production at HERA 

In the deep inelastic scattering regime, the heavy quark production at HERA is generally well described 
by NLO QCD. Massive scheme calculations combined with a conventional fragmentation function, fitted 
to data from e+e~ colliders, work well in describing c-meson production. Both the massive VFN and 
purely massive schemes describe the c and h production within the current uncertainties. Within their 
uncertainties, kt factorization calculations are also able to reproduce the data. 

In photoproduction, the descriptions by NLO QCD fall short in some regions. The predictions tend 
to undershoot the data in the region of lower transverse momenta and forward rapidity. Investigations 
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Figure 27: The first observations of liard diffraction by ZEUS (left) and HI (riglit). Left: 
Distribution in DIS events of rjmax, tfie maximum pseudorapidity of a cluster of energy larger 
than 400 MeV, in the ZEUS calorimeters, the range of which is indicated at the top. There 
is a clear excess seen at small rjmax of events in the data over a conventional DIS Monte Carlo 
simulation. Right: Similar observation in the HI DIS rjmax distribution, -E^m = 400 MeV. 



show that the discrepancy is largely in the resolved photoproduction regime and there are indications 
that missing higher orders is the cause. 

The data covered here are all based on about 120 pb~^ each for HI and ZEUS in the HERA I data 
taking period. There is about 400 pb~^ of data for each experiment which is yet to be analysed and 
published from the HERA II period. For ZEUS, the HERA II data set is taken with their newly installed 
microvertex detector. 

Due to space limitations, a large part of heavy quark studies at HERA could not be covered. These 
topics include studies of productions of c-hadrons other than the D* and the study of fragmentation 
functions |169j . the study of the transition region between DIS and photoproduction |170] as well as 
the study of the photoproduction region using the HI central silicon tracker \171\ . 



7 Diffraction 

7.1 Inclusive Diffraction 

Diffraction as a hard deep inelastic scattering process was discovered at HERA when ZEUS [172] and 
HI |173j observed an unexpected excess of DIS events characterised by the absence of activity in the 
forward direction which is usually occupied by hadrons emitted from the final state (Figure [27!) . In 
conventional QCD the hadronic final state is produced in the radiation from colour charges exchanged 
in the interaction. The observation of an excess of large rapidity gap (LRG) events, as the absence of 
hadronic activity at small angles is named, was readily understood as being due to a virtual photon- 
proton interaction without colour exchange. The basic process is sketched in Figure[28l An exchange 
of a colour singlet with vacuum quantum numbers is historically called Pomeron exchange. The early 
discovery of hard diffraction in DIS at HERA gave rise to the development of the chromodynamic theory 
of diffraction, based largely on the proof of collinear factorisation |174] of DDIS cross sections at fixed 
x/pand t. Here xjp is the fraction of proton momentum the diffractive exchange carries and t is the 
4-momentum transfer at the proton vertex. This allowed diffractive parton densities to be introduced 
with the constraint that the proton does not fragment during the violent ep collision. In fact the calculus 
of diffractive structure functions can be introduced quite analogously to the inclusive case [175] . 
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Figure 28: Schematic diagram of hard diffractive deep inelastic scattering. The scattered 
proton mostly stays intact but also may dissociate into a system of small mass My. The 
diffractive exchange carries a fraction xjp of proton's momentum. A parton of momentum 
fraction (3 = xjp/x interacts with the virtual photon which in DIS has a large virtual mass 
squared, Q^. The salient feature of these events is the absence of forward particle production, 
near the proton beam pipe, which is measured as a gap in polar angle, or equivalently rapidity, 
from the proton beam axis to the more centrally produced particles which stems from the 
struck parton in the diffractive exchange. Despite the violence of the interaction about a 
fraction of 10% (3%) is observed at of 10(1000) GeV^. 



Another way to analyse diffractive DIS is the colour dipole model (CDM) approach |176j in which, 
unlike in the parton model, the proton is studied at rest. Due to the smallness of x, the photon can 
be considered as disintegrating at a coherence length oc 1/x prior to interacting with the proton. This 
leads to a phenomenological prescription of diffractive scattering cross sections as the convolution of 
dipole wave functions and a dipole-proton cross section, which is to be determined from the data. 

The accuracy and scope of the diffractive measurements have improved dramatically over the past 
years. Much more accurate and complete data are now available based partially on the LRG iden- 
tification of diffractive DIS |177l 1178] . The agreement of all data is approaching an acceptable level, 
both for the LRG method data and also including more recent results applying the so-called Mx se- 
lection to diffractive scattering |179] . Inclusive diffractive measurements have been extended to high 

< 2000 GeV^. Diffractive DIS has also been observed in charged current scattering. The fraction of 
hard diffraction to inclusive DIS is about 10% at low ~ 20 GeV^ and diminishes, both in NC and 
in CC, to a few % at high ~ 2000 GeV^, where, however, the mean x is correspondingly larger by 
roughly two orders of magnitude at HERA. 

Both HI and ZEUS operated forward proton spectrometers in HERA I, tagging the proton down- 
stream, near the proton beam line with the "EPS" (HI) and "LPS" (ZEUS), at ~ +100 m distance 
from the interaction point. In the HERA II phase, HI installed in addition to the EPS a Very Eorward 
Proton Spectrometer, accessing an x/p range near to 0.01, while ZEUS had dismantled the LPS. The 
EPS and LPS detectors allow a direct comparison of the LRG based cross sections to the ones obtained 
with p tagging. Eor a diffractive mass range My < 1.6 GeV, HI measured the ratio of the LRG to the 
EPS cross sections as 1.23 ± 0.05 and to be constant within the kinematic range of the measurement, 
2 < < 30 GeV^, 0.01 < /3 < 0.7 and 0.001 < x/p < 0.05, where [3 = x/xjp is the momentum fraction 
of the diffractive exchange carried by its partons. The ZEUS result is 1.23 with a total error of 0.05 |180] 
for the ratio of the LRG to LPS data, in agreement with the HI measurement and also constant in the 
kinematic range of the measurement, 2.5 < < 40GeV2, 0.007 < /3 < 0.816 and 0.0002 < xjp < 0.02. 
Thus the dominant process in the gap data is indeed diffraction and further processes such as resonance 
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or neutron production do not alter the kinematic dependencies beyond the point to point uncertainty, 
of typically 5-10%, of the LRG/FPS or LRG/LPS ratio. 

It is satisfying to observe that the tagged differential diffractive cross section data are consistent 
between HI and ZEUS. The tagged data also allow the t dependence of the cross sections to be measured 
which can be parameterised as da/dt oc e~^*. A slope of 7 ± 0.3 is obtained by ZEUS |180j . HI has 
parameterised the slope as 6 = hip — 2aipliaxip corresponding to a trajectory a/p(t) = a/p(0) + ajpt. 
At small Xjpone obtains bjp = 6.0 with an error of about 2 and a small ajp = 0.02 consistent with 
zero within the errors |177] . A large a^p would have been an indication for shrinkage of the diffractive 
peak. 



7.2 Diffractive Parton Distributions 

The diffractive measurements based on the LRG method, which has the widest kinematic coverage, 
allowed QCD analyses to be performed in order to derive the quark and gluon distributions of the 
diffractive exchange. Such analyses are based on the diffractive cross section a^^^^ 

dxipdxdQ"^ xQ^ + r ( /i'' ); ( ) 

which is integrated over the ranges of t and My. Similarly to inclusive DIS, the reduced ep cross 
section depends on the diff: 
approximation according to 



section depends on the diffractive structure functions and F^*-^-* in the one-photon exchange 



For y not too close to unity, a^^^^ = F^^^^ holds to very good approximation. 

Recently a detailed QCD analysis was published by HI |177j . Since data are used for a range of 
x/p values, it is required to factorise the x/p dependence out which in |177j is done with an ansatz 
inspired by Regge theory. The QCD fit, similar to truly inclusive scattering as described above, then 
determines parton distributions, here a singlet quark distribution I]{x,Q'^), assuming u = d = u = 
d = s = s, and a gluon distribution xg{x,Q^), where x = Pxjp. The data, as is shown for example 
in Figure [29] (left) for a medium x/p value, are well described in their x and dependence by this 
approach. In the analysis, rather stable results were obtained for both distributions, with the exception 
of the behaviour of the gluon distribution at large momentum fractions which is rather uncertain. This 
can be understood from the consideration of the In derivative of the diffractive reduced cross section 
(see Figur e[29l (right)). As the two QCD fit components reveal, at large j3 there is not much sensitivity 
to gluons Thus HI decided to publish two fits, with very similar but different xg at large /3. 
A distinction between these two possibilities was made possible with diffractive dijet data |181j . The 
diffractive dijet production leads to improved sensitivity to the gluon distribution in the diffractive 
exchange (see Figure lSUl) . Thus the gluon distribution uncertainties can be reduced significantly, similar 
to the influence of dijet DIS data on the determination of the gluon distribution in the proton. 

There are many further similarities between diffractive and inclusive DIS: The logarithmic derivatives 
of (7^ and ar, with respect to In Q^, are measured to be the same for (3 < 0.6. Thus low (3 diffraction and 
low X inclusive scattering seem to reflect a common origin, the dynamics of the QCD vacuum. It is not 
surprising, then, that a number of recent observations in diffractive DIS resemble those in inclusive DIS; 



^This result is reminiscent of inclusive DIS: at large x the scaling violations are driven by only and no direct 
sensitivity to the gluon distribution appears. This is the main reason why older fixed target DIS experiments have not 
determined xg well. It required a high energy ep collider to measure the gluon distribution in the proton. The same holds 
for polarised ep scattering, the determination of the polarised gluon distribution G cannot reliably be obtained from lower 
energy data as these are essentially sensitive to chromodynamic bremsstrahlung only. 
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Figure 29: Left: Diffractive reduced cross section at xjp =0.03 as a function of for different 
j3 = x/xip. The data are HI data from HERA I. The curves are different NLO QCD fits, 
dashed: fits A and B (see text) to the inclusive data only, solid: fit to inclusive and diffractive 
dijet data. Right: Derivative of the reduced diffractive cross section at xip = 0.01 and its 
decomposition into the gluon and quark parts in the HI NLO QCD analysis. 

such observations include the rise towards low /3 or x of the In derivatives and the fraction of charm 
of about 20%, for away from thresholds, in both the diffractive cross section and in the inclusive 
cross section. Diffraction has evolved from the first surprise of its presence at HERA to a rather exact 
and interesting testing ground of perturbative QCD. Further data from both HERA collaborations and 
studies of universality in comparing DIS ep and photoproduction predictions |183j . as well as of HERA 
data with Tevatron and later LHC data, will certainly contribute to this field in the future. Recently, 
diffraction is being considered as a means of studying possible SM and SUSY Higgs particle production 
at the LHC, in a probably rare but possibly clean double diffractive scattering process pp pHp [184j . 

7.3 Colour Dipole Model Description of Diffraction 

Diffractive parton densities lead to predictions in the same way conventional parton densities do. How- 
ever, the xip dependence, which is some measure of the colour singlet formation within the proton, is 
factored out and remains a phenomenological assumption made prior to the subsequent QCD analysis. 

An alternative way to study the phenomenon of diffractive DIS is through the application of the 
Colour Dipole Model (CDM). The cross-section is here given by the convolution of the quark-antiquark 
colour dipole wave-function and the dipole-proton cross-section. This formalism connects the inclusive 
low x cross-section, discussed in Section 3.4.3, to the DDIS cross-section in the same way forward elastic 
hadron scattering is related to the total cross-section. 

The dipole-proton cross-section may be describable in terms of the DGLAP, or the BFKL formalisms 
in their respectively applicable kinematic regions. This has not yet been fully demonstrated. In the 
very low x or low regions, gluon recombination and/or non-perturbative effects should influence the 
cross-section. 

Figure [3T] shows an example of a fit to DDIS cross-section data [185 J using a CDM-type model, 
BEKW [186j . The description achieved can be very good, and some models describe a number of cross- 
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Figure 30: Diffractive partons as determined in various NLO QCD fits: dotted: inclusive HI 
fit (A); dashed: inclusive HI fit (B); solid: HI fit to inclusive and dijet data; dashed-dotted: 
from [182| . Here z is the momentum fraction of the diffractive exchange carried by the 
parton. 



sections, DDIS, inclusive DIS, Heavy Quark production, simultaneously. While the phenomeno logical 
success is remarkable, it is fair to say that at this stage, there is no theoretical consensus in the 
interpretation of these models. The reader is referred to recent articles and reviews pUl I176j for further 
discussions. 

7.4 Elastic Vector Meson production 

The description of exclusive electroproduction of vector mesons, ep eVp, is closely related to the CDM 
model described above. The virtual photon emitted from the electron fluctuates into a quark-antiquark 
dipole, which then scatters elastically off the proton. Long after the interaction, the qq pair forms a 
vector meson. At sufficiently high energies, the three steps, the qq formation, the dipole scattering 
and the formation of the vector meson are well separated in time. The cross-section for the process, 
then, is factorised into the qq coupling to the photon, the dipole-proton scattering amplitude, and the 
final-state formation [1871 1188] . The first and the last steps involve the electromagnetic coupling of the 
qq pair and the wave function of the vector meson, respectively. In case of light mesons, parton-hadron 
duality is sometimes invoked for the last step. 

If there is a sufficiently large scale involved, the dipole-proton cross-section is calculable in terms 
of two-gluon exchange. In principle this cross-section depends on the generalised parton distribution 
function (GPD) |189j of the gluon. GPDs contain information on the correlation of partons within the 
proton. At low x, i.e. where leading Inl/x approximation is valid (and if t is small), the generalised 
gluon distribution can be approximated by the usual gluon distribution. The cross-section then should 
rise steeply in 1/x (or equivalently in W, the virtual photon-proton cms energy, which at low x is given 
by ^s/x.) for fixed refiecting the steep rise of the gluon density at low x. 

For those cases where the interaction is soft, the cross-section should rise slowly with W, similarly 
to the well-known Regge energy behaviour of hadron-hadron total cross-sections. 

Figure [521 left) show the measurements of elastic VM photoproduction |191j from HERA as well as 
the fixed-target experiments. The data are fitted to the form to quantify the rate of the rise of 
the cross-section with energy. The 6 increases as the mass of the VM increases. For the J/\E' one finds 
a value of 5 ~ 0.8 which is consistent with a gluon density varying as xg{x) oc close to typical 



52 



BEKW(mod): — Total ■■■■ (qq)^ 
P=0.005 P=0.025 P=0.125 P=0.4 



(qq)L 

P=0.7 



(qqg)T 

P=0.9 



0.05 




10 1 

(GeV^) 



Figure 31: DDIS measurements by ZEUS |179j . The lines are fits to a CDM-type model 
BEKW [T86] . 

exponents extracted from the QCD fits to DIS inclusive data at small Qq ~ Mj^^. 

Figure [32^ right) shows the quantity 6 extracted from many different measurements |190l I191j of 
VM production at HERA, now as a function of + My, where My is the mass of the vector meson 
being produced. While the uncertainties are still large, at low values of + My, S ~ 0.2 which is the 
expected value for Regge energy behaviour while at higher + My the exponent 6 rises to around 0.8. 
The data are consistent with + My being the appropriate hard scale. 

7.5 Deeply Virtual Compton Scattering 

The process ep —>■ e^p at the HERA collider has contributions from both the Bethe-Heitler (BH) and 
the deeply virtual compton scattering (DVCS) processes (Figure 1551) . In the kinematic region of the 
HERA collider, the interference term between DVCS and BH processes is small, and the latter, which 
is precisely calculable, can be subtracted from the measured cross-section in order to obtain the DVCS 
cross-section. 
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Figure 32: (left) Elastic VM photoproduction cross-section as a function of W for a variety 
of species of VMs for HERA and fixed-target experiments. The data are fitted to the form 
W^, and the results are shown in the plot, (right) The extracted values of 6 from HERA 
elastic VM production and DVCS data as a function of + My^j. 



The DVCS process is rather similar to the elastic VM production process, but has the advantage 
of having a simple final state with no need to rely on an understanding of meson wave functions. 
Figure [31] shows the measured DVCS cross-section |192l I193j as a function of W . The rise with W is 
relatively steep (5 ~ 0.75, see Figure [321) demonstrating the relative "hardness" of this cross-section. 
The theoretical predictions of Freund and McDermott [194] use GPDs based on normal PDFs (MRST 
and CTEQ) [T95| . 

In a class of CDM models called the saturation model |197j . the dipole-proton cross-section is 
proportional to the square of the dipole separation at small r oc (l/Q). At sufficiently large r, 
the dipole-proton cross-section saturates (becomes constant as a function of r) to preserve unitarity. 
The transition of a process into the saturation regime is characterised by the saturation scale Qs{x) = 




Figure 33: Diagrams illustrating the DVCS (left) and the Bethe-Heitler (center and right) 
processes. 
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Figure 34: The DVCS cross-section as a function of W for Q"^ = 8 GeV^ and |t| < 1 GeV^. 
The inner error bars are statistical and outer, statistical and systematic uncertainties taken 
in quadrature. 



Qo{xo/x)~''^^'^ , where Qq, Xq and A are parameters. One prediction of such models is that certain DIS 
cross-sections can be expressed as a function of the single variable r = Q'^ /Ql{x). Figure [35] (left) is 
the total ep DIS cross-section data at low-s (< 0.01) as a function of r |196j . In this case Q'^ = 1 
GeV^, xq = 3.04-10-^ and A = 0.288, as was determined in [M]. Figure [35] (right) show the HI DVCS 
data jl93j using slightly different parameters, Ql = 1 GeV^, Xq = 2.7-10"^ and A = 0.25. It has been 
shown that such a scaling holds relatively well also for elastic VM production and DDIS cross-sections. 



8 Electroweak Measurements 
8.1 Charged Current Cross-Section 

Rewriting Equation [231 the charged current DIS cross-section (Pa^'^ /dxdQ^, for the reaction ep — *• 
i'{u)X, can be expressed as: 

dxdQ^ 27rx ^M2, + Q2J ^r-,cc, l^yj 

where cr^cc contains three charged current proton structure functions depending on the lepton beam 
charge. If the Fermi constant Gp, is equated with that measured in muon decay, then the residual 
correction to the above equation from higher order electro-weak corrections amounts to only a few 
parts per mille [198J. 

If the structure functions are known, then the measurement in the region -C is a measurement 
of Gp independently of Q^. In the region ^ M^r, the dependence of the cross-section is a 
measurement of the mass of the particle exchanged in the t channel, namely the W. 

Figure [36] shows the charged current cross-section as a function of at HERA for all data |199] 
with unpolarised initial leptons. The effect of the propagater mass can be seen as increases towards 
M^. The difference between e~^p and e~p charged current cross-sections arises mostly from the fact 
that in the former primarily the (i-valence quark is probed whereas in the latter the w-valence is probed. 
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Figure 35: Left) Total DIS cross-section (written as the 7*p cross section) as a function of the 
saturation scahng variable r. The figure is taken from |196j . Right) The DVCS cross-section 
from the HI collaboration function of r. 



The neutral current cross-section [1991 I200j is also compared to that of the charged current, in 
Figure [361 It is observed that the two cross-sections indeed become about the same at the scale of 

^ 10^ GeV^ giving an explicit demonstration of the electro-weak unification as had been anticipated. 

A unique check of the standard model can be made by using the precisely measured muon decay con- 
stant (1.1639- 10~^ GeV~^ [32]) and determining the propagator mass M]y that appears in Equation [^ 
A fit, with as a free parameter gives, for example: 

Mw = 78.9 ± 2.0(stat.) ± 1.8(syst.)l?;0(PDF)GeV, (50) 

for the ZEUS determination using e^p data from HERA I only. The HI collaboration determined: 

Mw = 80.9 ± 3.3(stat.) ± 1.7(syst.) ± 3.7(theo.)GeV, (51) 

using 35.6 pb~^ of e"^p data from HERA I. These uncertainties are much larger than the direct measure- 
ments of Mw and the measurements in the time- like region [12]. However, the HERA measurements 
demonstrate a consistency of the space-like measurements with the Standard Model. 

The purely weak nature of the charged current interaction means that the cross-section for polarised 
initial leptons is directly related to the unpolarised cross-sections as: 

c^'<?(e^p) _ ^ p^ d^(^unpoiie^P) ^52) 
dxdQ"^ dxdQ"^ 

where the longitudinal polarisation of the lepton beam is defined as 

„ Nr-Nl 



Nr + Nl 



(53) 



where Nr and Nl are the numbers of right- and left-handed leptons in the beam, respectively. 

With the advent of longitudinally polarised lepton beams in collider mode, in HERAH, it became 
possible to directly measure the polarisation dependence of charged current DIS. The result [1991 1201[ 
I202j is shown in Figure [371 
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Figure 36: Left: The charged current cross-section as a function of measured at HERA 
I. Both the e^-p and e~p cross-sections are shown. Right: the charged current cross-sections 
are compared with the neutral current cross-sections measured at HERA. 



Within the standard model, Gp is a function of the fine structure constant, a, the masses of the 
gauge bosons M^y and Mz- Through radiative corrections, Gp also depends on the top mass, Mf, as 
well as on the Higgs mass, M/^, though only logarithmically. Rather than taking Gp from the muon 
decay constant, it is possible to take values of a, Mz and Mt from other experiments, and determine 
My^ at HERA. Given the uncertainty of Mt from the Tevatron, which is at the level of 2 GeV [203], it 
is estimated that the uncertainty of Mvi/ in the order of 50 MeV will be possible using the full statistics 
from HERA [204j . This uncertainty, which is comparable to those of current direct measurements |42j . 
is not enough to predict the Higgs mass with any certainty; however it is interesting to see if the trend 
at the Tevatron for M^ and Mj measurements to predict a low mass for the standard model Higgs 
boson is confirmed by this result. 



8.2 Neutral Current Cross Section 

A polarisation asymmetry measurement of the NC scattering, for helicities Pr and Pl, according to 
EquationsfTSl and [191 determines a combination of F^^ and xF^^ 

-TR — -fL y+ 

neglecting the pure Z exchange terms, which are small at HERA. The second term is a small correction 
since the vector coupling is small, as is the factor YL = 1 — (1 — y)^, in most of the kinematic range 
at HERA. The product cte-F^^ is proportional to combinations a^Vq and is thus a direct measure of 
parity violation at very small distances, ~ 10~^^ m, as they are probed with electroweak cross section 
measurements at HERA. 
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Figure 37: The total cross-section for e^p CC DIS with > 400 GeV^ as a function of 
the longitudinal polarisation of the lepton beam. The lines show the prediction of the SM 
evaluated using CTEQ and MRST PDFs. 



The polarisation NC cross-section asymmetry 



to a very good approximation measures the structure function ratio 

~ ^Kzae^. (56) 



Thus A"*" is expected to be positive and about equal to —A . At large x these asymmetries measure 
the d/u ratio of the valence quark distributions according to 

The preliminary measurements of from HERA, Figure [3H1 based on the combined data |201j from 
HI and ZEUS show a significant polarisation effect and the asymmetries A^ to be of opposite sign 
as predicted. The lines, which describe the data well, are the predictions of the standard model as 
obtained from the HI and ZEUS QCD fits. 

It is clear from Equations [18] to EU that neutral current cross-sections at HERA are sensitive to the 
vector and axial- vector couplings of the quarks. The sensitivity of the cross-sections to through xFs is 
enhanced by polarisation. Unpolarised cross-sections are sensitive to Vg mainly via F2 whereas polarised 
cross-sections have further sensitivity through XF3. Figure [39] shows the published measurement of Og 
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Figure 38: Measurements of the polarisation asymmetries by the HI and ZEUS collab- 
orations. The error bars denote the total uncertainty. The curves describe the standard 
model predictions using the PDF sets obtained from unpolarised data. 



and Vq by the HI collaboration compared to determinations from LEP and Tevatron. The figure also 
shows the preliminary determinations |2U5j of ZEUS using the polarised data which shows the expected 
improvement in precision. 

From Equations [18] and [191 it can be seen that the difference between the unpolarised e~^p and 
e~p cross-sections can be defined to be xFs. Terms proportional to Ve, which is small (0.036), can be 
neglected so that, to leading order in pQCD, one has 

xFf ~ ±aeKzxF;^^ (58) 

to a good approximation, with 

xFf = 2x[eA{U -U) + e^D - D)] (59) 

and U = u + c and D = d + s for four flavours. Assuming that there are no anomalous differences 
between sea quark and anti-quark distributions, xF^^ depends on the valence quark distributions of 
the proton only, and on the axial vector couplings of the u and d quarks. 

The preliminary measurement [201] by the ZEUS and HI collaborations of this quantity as a function 
of X at of 1500 GeV^ is shown in Figure HOI The measurements are well described by the standard 
model predictions using PDFs as obtained from the HERA data. 

Substituting the standard model values for the couplings, one finds 

xFf = |(2m„ + 4) (60) 

where and d^ are the valence up and down quarks distributions. To leading order, the above equation 
leads to a sum rule [206] 

f xFf— = \ [\2u, + d,)dx = I. (61) 
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Figure 39: The 68% confidence level contours of the electroweak parameters vs Vu and 
vs Vd from the fits to HI (polarised only) and preliminary ZEUS (polarised and unpolarised 
data) compared to the results from LEP and Tevatron. 



In the range of the data, the integral of F^'^ is measured to be 

fO.65 
/0.02 



/■0.65 

/ Ffdx = 1.21 ± Om{stat) ± 0.08{syst) (62) 

J0.02 



which is in agreement with predictions of the SM for this x range, using, for example, HI and ZEUS- Jets 
PDFs. 

9 Searches 

Searches for "Beyond the Standard Model" (BSM) phenomena at HERA is a large topic, reviewed at 
some length in (20 7j relatively recently. Here only a few selected topics are covered based on updated 
results. 

Unlike LEP and the Tevatron, HERA does not have a particle and its anti-particle in the colliding 
beams. Therefore, pair production cross-sections of particles are rather small but new particles may 
favourably be produced singly. As a result, the limits obtained at HERA are generally presented as a 
function of the couplings of the new particle as well as its mass. 



9.1 Leptoquarks 

In general, measurements for BSM searches may be interpreted, or analysed, in the context of a par- 
ticular model ("model-dependent analysis"). On the other hand, the data, very often the very same 
data used in the mo del- dependent analysis, can be interpreted simply as possible deviations from the 
Standard Model. A first example of the former is the search for Leptoquarks. 

Leptoquarks (LQs) [208], colour triplet bosons carrying both baryon and lepton numbers are pos- 
tulated in many extensions of the Standard Model. An example are squarks in R-parity violating 
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Figure 40: The measurement of the structure function xF^^ by the HI and ZEUS collab- 
orations. The inner error bars denote the statistical uncertainty while the full error bars 
comprise the statistical and systematic uncertainty added in quadrature. The data are from 
HERA I but include also part of the polarised HERA II data, in which, however, effectively 
polarisation effects average out leading to a small correction only, see [201] . The curves 
describe the Standard Model predictions as obtained in the HI and ZEUS NLO QCD fits. 



supersymmetry |209j . In ep collisions at HERA, LQ states may be produced directly though e-quark 
fusion, or in the s-channel, with subsequent decay into eq or z/g yielding resonant peaks at the LQ 
mass. LQs may also be exchanged in the w-channel, in which case characteristic deviations from the 
SM cross-sections are expected at high Q^. 

Both HI and ZEUS have published the results of searches |21UII211] for LQs from the HERA I period 
using 120-130 pb^^ of data each. In the HI search for LQs, the reconstructed CC and NC events are 
searched for deviations from the SM. In the s-channel production of LQs, an enhancement is expected 
at the X corresponding to the LQ mass of Mlq = ^yxseq- In the ZEUS search, the invariant mass of 
the final-state e-jet and z/-jet system is reconstructed directly; in this case Mlq = Me{u)-jet- The two 
methods are rather similar in the large region which corresponds to large LQ masses. 

Neither search finds evidence for LQs and limits are set. At HERA, it is usual to set limits based on 
the phenomenological model proposed by Buchmiiller, Riickl and Wyler (BRW) |208j which describes 
fourteen types of LQs. It is out of the scope of this review to give the limits for each of these. As an 
example. Figure W2\ shows the limits obtained by the experiments for one of the 14 BRW LQs compared 
to the limits set at LEP and the Tevatron. At LEP, the searches are sensitive mostly to the coupling, 
A, as is the case at HERA above the kinematic limit of M = ^/s since both rely on the virtual exchange 
of the LQs. The Tevatron searches, on the other hand, are mainly sensitive to the mass since the LQs 
would be pair produced via quark-anti-quark annihilation. 

There are many ways to state the limits, and the reader is referred to the original papers for a 
fuller discussion. Also limits for LQ couplings to higher generation quarks have been set in searches for 
lepton-fiavour violation by both ZEUS and HI collaborations |216l 1217] . 



Hl+ZEUS Combined (prel.) 

Q^=1500 GeV^ 
— HI 2000 PDF 
---- ZEUS- JETS PDF 
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Figure 41: Diagrams for s-channel leptoquark production and w-channel LQ exchange. Here 
e stands for both electron and positron, and v for both neutrino and anti-neutrino. 



9.2 Excesses Beyond the SM 

Rather than doing a mo del- dependent search for phenomena beyond the Standard Model, one may study 
processes with a low cross-section in the SM in order to observe possible deviations from expectations. 

In 2002, the HI collaboration reported an excess of events with high energy isolated electrons or muon 
accompanied by missing transverse momentum |218j . In approximately 100 pb~^ of data, HI observed 

10 events with transverse energy of the hadronics system, P;^, greater than 25 GeV. Such events are 
expected in the SM from W boson production with subsequent leptonic decay. The prediction from the 
SM and background from misidentified CC and NC events was 2.9±0.5. In a similar search, the ZEUS 
collaboration found no excess over the SM |219j . 

Both collaborations have recently made preliminary updates |220] including the data from HERA 

11 which amounts to a five fold increase in integrated luminosity. Figure US] shows the distributions of 
these events as a function of as reported by the HI collaboration. 

The results from both collaborations in the region of excess observed by HI originally is summarised 
in Table HI Overall, the data show a good agreement with SM expectations. In detail, there is an 
excess observed by HI in the e^j) sample for both electrons and muons. The excess is at the level of 3 
a if these two HI samples are summed. The ZEUS data in the same channel, on the other hand, is in 
good agreement with the SM. In a combined analysis [221j . where a common phase space was defined, 
29 isolated electrons and muons were observed at > 25 GeV, in 0.97 fb'^ of data, where 25.3±3.2 
events are expected from SM. Both HI and ZEUS collaborations have also looked in the r channel, and 
have found no significant excess beyond the SM [222[ 1223] . 

Having done a model-independent search, one can use the data to set limits on particular models. 
One interpretation of this type of excess would be an anomalous coupling of the u-quark to the t-quark. 
Such a coupling in the SM is negligibly small. Indeed the same type of analysis has been used to set 
limits on a single top production at HERA, and thus on the anomalous twy and tuZ couplings [216U224j . 

The HI collaboration had observed events with high Pt di- or tri-electrons in excess of those expected 
by the SM via higher order electroweak events [225] : there were three di-electron events and three tri- 
electron events, with the mass of a pair of electrons in excess of 100 GeV, in a data sample corresponding 
to the integrated luminosity of 115 pb^^ where the expectations were 0.30±0.04 and 0.23±0.04 events 
respectively. 

In recent updates, both HI and ZEUS collaborations have released preliminary results on multi- 
electron production using over 450 pb~^ each |226] . Neither collaboration finds a significant excess 
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Figure 42: Exclusion limits at 95% confidence level on the coupling A as a function fot lepto- 
quark mass for Sq^l in the framework of the BRW model from HI and ZEUS collaborations. 
The limits from L3, OPAL and DO collaborations [2l2l [2131 MM are also shown. The DO 
limit in the ZEUS plot is an older one |215] current at the time of the publication of the 
ZEUS result. 

beyond the SM in the larger data sample. HI collaboration searched for multi-lepton events, where 
leptons are either electrons or muons, and found four events with a scalar sum of lepton transverse 
momenta greater than 100 GeV, whereas 1.9ib0.4 events are expected from SM. ZEUS searched for di- 
and tri-electrons with the mass of a pair of electrons in excess of 100 GeV, and found 6 events where 
5.4±1.1 events are expected from SM. In a combined analysis |22 7] , done in a common phase space, 
HI and ZEUS together observed 6 di- and tri-electron events with a scalar sum of electron transverse 
momenta greater than 100 GeV, in 0.94 fb^^ of data, where 3.0±0.3 events are expected from SM. 

10 Summary 

HERA was the first electron-proton collider ever built. It allowed the investigation of physics of deep in- 
elastic lepton hadron scattering to be extended by two orders of magnitude into a new kinematic domain 
in terms of four-momentum transfer squared, Q^, and Bjorken x. Unlike previous fixed target experi- 
ments, which scattered either charged leptons (electron and muons) or neutrinos off nucleon targets, the 
very large region at HERA, at which weak interactions become manifest, allowed the simultaneous 
investigation of neutral current {ep — eX) and charged current {ep vX) interactions. Furthermore 
almost real photon-proton interactions - the photoproduction processes - could be investigated in which 
particles or jets at high masses are produced. 

The accelerator, comprising separate and p sources, pre-accelerators and rings, was a new chal- 
lenge for machine physics. The challenges included the problem of polarised ep collisions, which was 
mastered successfully and will provide the a knowledge base for future ep machines such as that being 
considered in conjunction with the LHC p beam. The collider detectors, HI and ZEUS, were developed 
independently and used very different detector techniques; for example, calorimeter techniques were 
based on liquid argon and uranium at HI and ZEUS, respectively. The physics results obtained turned 
out to be of similar quality and coverage and in good agreement. Occasional discrepancies between 
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Table 4: The summary of HERA preliminary results for events with isolated electrons or 
muons and missing transverse momentum with > 25 GeV. The number of observed 
events is compared to the SM prediction. The signal component of the SM expectation, 
dominated by W production, is given as a percentage in the parenthesis. The selections for 
HI and ZEUS results are similar but not identical. See |220j for details. 



Isolated e candidates 


HI Prelim. 


ZEUS Prelim. 


e'p (HI 


: 184 pb-i, ZEUS: 204 pb"!) 


3/3.8 ±0.6 (61%) 


5/3.8 ±0.6 (55%) 


e+p (HI 


: 294 ph-\ ZEUS: 228 pb-^) 


11/4.7 ±0.9 (75%) 


1/3.2 ±0.4 (75%) 


e^p (HI 


: 478 ph-\ ZEUS: 432 pb^^) 


14/8.5 ± 1.5 (68%) 


6/7.0 ± 0.7 (64%) 




Isolated /i candidates 


HI Prelim. 


ZEUS Prelim. 


e^p (HI: 


184 pb-\ ZEUS: 204 pb^^) 


0/3.1 ±0.5 (74%) 


2/2.2 ±0.3 (86%) 


e+p (HI: 


294 ph-\ ZEUS: 228 pb-^) 


10/4.2 ±0.7 (85%) 


3/3.1 ±0.5 (80%) 


e^p (HI: 


478 pb-\ ZEUS: 432 pb-^) 


10/7.3 ± 1.2 (79%) 


5/5.3 ±0.6 (82%) 



the data or methods used led to constructive and sometimes challenging interactions between the two 
collaborations. 

The intra-Hl-ZEUS collaboration is currently being moved to a new level, in which HI and ZEUS 
will combine their final data sets where appropriate. An example is the combination and common QCD 
analysis of the inclusive NC and CC cross section data. Such an analysis can be used to, in essence, 
cross-calibrate the results, taking advantage of the different detector techniques and different kinematic 
reconstruction methods. The combination is expected to lead to results of a new quality, exceeding 
what is to be expected in accuracy from a simple statistical average of the data. HERA will, in this 
way, provide a base for predictions of various processes to be studied at the LHC; and QCD will have 
to be used to extrapolate HERA measurements to the LHC kinematic regime. While the final results 
of HERA are being obtained, the first observations at the LHC are expected to be made. This will lead 
to a fruitful link between ep and pp physics, between HERA and the LHC. 

The main results of HERA as had been presented in this overview may be summarised as follows: 

• HERA opened the field of low x physics having reached values of x = Q'^/sy of the order of 
10~^ at larger than a few GeV^, owing to the unprecedented large energy s. In this newly 
accessed range, parton densities were discovered to be rising towards lower x at fixed which 
is the consequence of gluon dominance at low x. This is in contrast to valence quark dominance 
at high X. Thus, a new phase of matter was discovered in which the densities are high but the 
coupling constant is small compared to unity. This discovery has immediate consequences for 
the understanding of quark-gluon dynamics in QCD as well as for high density states which are 
being investigated, for example, in nucleus-nucleus scattering or in super-high energy neutrino 
interactions. 

• A further discovery of HERA was the observation of hard diffractive scattering comprising a sig- 
nificant fraction (order 10%) of the cross-section. This, attributed to the exchange of vacuum 
quantum numbers, is a process subject to perturbative QCD calculations due to the large mo- 
mentum transfer squared involved. These colourless exchange reactions may lead to a firmer 
understanding of confinement, and, at the LHC, to the measurement of quantum numbers of 
the Higgs boson in a supersymmetric extension of the theory, in a rather clean experimental 
environment. 

• As a machine of cleanest resolution of the structure of matter, HERA has set a new limit of 
7 • 10~^^ m to a possible substructure of quarks, about 10'^ times smaller than the proton radius. 
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Figure 43: The distribution of transverse hadronic energy, P^, for events with isolated 
leptons as measured by the HI collaboration. The e^p (left) and e~p (right) data are 
compared to the SM expectations (histogram labeled "All SM"). The component due to 
true SM signal, dominated by W production, is given by the histogram labeled "Signal". 
^data is the total number of data events observed, ^sm is the total SM expectation. The 
total uncertainty of the SM expectation is given by the shaded band. 



While one could have hoped to find a preonic substructure, this result constitutes an important 
milestone and signals that higher energy experiments will have to be performed in order to look 
even deeper into the structure of matter. 

• The results of HI and ZEUS have led to an unexpectedly rich harvest of QCD related results which 
are the basis of developing a further understanding of the theory of strong interactions. This refers, 
for example, to the findings on multi-jet production in DIS, the measurement of the strong coupling 
constant and the determination of a rather complete set of parton distributions, the analysis of 
heavy quarks as dynamically produced in boson-gluon fusion, and the first measurements on 
deeply virtual Compton scattering at low x which access parton correlations. 

• The recent results of HERA have a precision which is becoming competetive for certain measure- 
ments of electroweak theory parameters at highest energies; an example is the determination of 
the light quark weak neutral current couplings. 

There have been many searches performed for signals of new physics, new particles, new symmetries, eg 
resonances and QCD phenomena. At the present stage of the data analyses, a few years prior to the final 
results, there have been no signals observed of SUSY phenomena, R-parity violating single production 
of SUSY particles, leptoquarks or high density instanton states. The limits set on new particles are 
competitive and complementary to the limits set in the crossed amplitude reactions at similar fermion 
energies, as have been studied in e"*"e~ reactions at LEP and in pp reactions at the Tevatron. If there 
is physics beyond the Standard Model, near the accelerator energy frontier, it is thus most likely to 
appear at TeV energies rather than at energies of a few hundred GeV, which characterise the Fermi 
scale as has been thoroughly studied in the last two decades of particle physics. 
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The operation of HERA ended in 2007. The analysis of data taken by the colhder experiments ZEUS 
and HI, at the time of this pubhcation, is expected to still take a few years. More precise results and 
also new results will still emerge. The remaining puzzles, such as that on the existence of penta-quark 
states, is expected to be resolved. The physics of deep inelastic ep scattering has proven to be a vital 
part of particle physics, and it may possibly be continued at even higher energies than HERA. 

Acknowledgements 

The overview presented here is a brief summary of the work of perhaps a thousand technicians, engineers 
and physicists for more than two decades. The authors had the great privilege to be part of a most 
exciting period of particle physics which in a few years will come to an end. It is our great pleasure 
to thank many colleagues, too numerous to be named here, for the possibility to share their wisdom 
and expertise and for most fruitful years of collaboration. As spokesmen of HI and ZEUS during the 
times where HERA's upgrade was in some difficulty, it is particularly pleasant to acknowledge the 
success of the HERA crew and to thank them for their impressive efforts without which the collider 
experiments could not have been operated so well. We acknowledge the strong support of the HI and 
ZEUS experiments by the theoretical particle physics community. DESY had been a generous host to us 
and many of our collaborators from many countries. It is not forgotten, that the wisdom and dedication 
of Volker Soergel, Bjoern Wiik and others laid the foundations so that the following generation (to which 
we belong) was able to work on a most exciting machine and on some of the most fundamental problems 
in elementary particle physics. 

We would hke to thank M. Derrick, S. Magill and P. Nadolsky for careful reading of parts of this 
manuscript. Any remaining inconsistencies or errors are, of course, our own. We would like to thank the 
HI and ZEUS collaborations for permission to show some of their results, presented so far at conferences 
only, also in this review. 

RY's work at Argonne National Laboratory was supported by the U.S. Department of Energy, Office 
of Science, under contract DE-AC02-06CH11357. 

References 

[1] B. Wiik and C. Llewellyn Smith, "Physics With Large Electron- Proton CoUiding Rings," preprint 
DESY 77/38 (1977) (unpublished) 

[2] HERA and the LHC, Workshop Proceedings, CERN-2005-014, DESY-PROC-2005-014, eds. 
A. DeRoeck and H. Jung, see also http://www.desy.de/~heralhc/ for further studies. 

[3] HI Collaboration, sec http://www.hl-desy.de 

ZEUS Collaboration, see http://www.zeus-desy.de 

[4] "Deep Inelastic Scattering and QCD", DIS Workshop Series. Since 1994 this annual series of 
workshops provides a most complete overview on HERA physics. For recent workshops see: 
DIS07, Munich, Germany, April 2007, proceedings to be pubhshed, http://dx.doi.org/, ed. 

G. Grindhammcr and K. Sachs 

DIS06, Tsukuba, Japan, April 2006, proceedings, eds. M. Kuze, K. Nagano and K. Tokushuku 

(2006) 

[5] H. Abramowicz and A. Caldwell, Rev. Mod. Phys. 71 (1999) 1275 

[6] R. Devenish and A. Cooper-Sarkar, "Deep Inelastic Scattering" , Oxford Univ. Press (2004) 

[7] E. D. Bloom et al., Phys. Rev. Lett. 23 (1969) 930 
M. Breidenbach et al., Phys. Rev. Lett. 23 (1969) 935 



66 



[8] BCDMS Collaboration, A.C. Benvenuti et al., Phys. Lett. B 223 (1989) 485 

[9] EMC Collaboration, J. J. Aubert et al., Phys. Lett. B 123 (1983) 275 and Nucl. Phys. B 259 (1985) 
189 

[10] A. Febel, H. Gerke, M. Tigner, H. Wiedemann and B. H. Wiik, IEEE Trans. Nucl. Sci. 20 (1973) 
782 and DESYF35,H,Sl-73/l 

[11] PEP Conceptual Design Report, LBL-4288/SLAC-189 (1976) 

Particle Physics with CoUiding Proton-Electron- Positron Beams, LBL-750/SLAC-146 (1972), (un- 
published) 

[12] T. Kamae and J. Shimizu, eds., Tristan ep Working Group Report, University of Tokyo, UTPN-165 
(1980), (unpublished) 

[13] B. Wiik et al, PROPER, DESY 1974, (unpublished) 

[14] J. R. Ellis, B. H. Wiik and K. Hiibner, "CHEEP: An ep Facility in the SPS," CERN- 78-02 (1978), 
(unpublished) 

[15] W.R. Frisken et al., "CHEER Canadian High Energy Electron Ring Feasibility Study", Fermilab 
Proposal 703 (1980), (unpublished) 

[16] HERA Technical Proposal, DESY-HERA-81-10 (1981), (unpubhshed) 

[17] A. A. Sokolov and I. M. Ternov, Sov. Phys. Dokl. 8 (1964) 1203 

[18] HI Collaboration, Technical Proposal, Hamburg 1986, (unpubhshed) 

[19] ZEUS Collaboration, Technical Proposal, Hamburg 1986, (unpublished) 

[20] HERMES Collaboration, see for example the gi measurement: A. Airapetian et al., Phys. Rev. D 
75 (2007) 012007 

[21] HERA-B Collaboration, see for example the T cross section measurement: I. Abt et al., Phys. Lett. 
B 638 (2006) 13 

[22] V. Andrecv et al., HI Background Working Group, "Technical Report on the Beam Induced Back- 
grounds in the HI Detector", HI Note 606/10(2002) (unpublished). Available upon request to 
HI 

[23] D. Bailey et al., ZEUS Background Working Group, "Study of Beam-Induced Backgrounds in the 
ZEUS detector from 2002 HERA running", ZEUS note 02-018, addenda 02-20 and 02-27 (2002) 
(unpublished). Available upon request from ZEUS 

[24] T. Alexopoulos et al, "eD Scattering with HI, A Letter of Intent", DESY 03-194 
H.Abramowicz et al, "A New Experiment for HERA", MPP-2003-62 
M. Klein and T. Greenshaw, J. Phys. G Nucl PaH. Phys. 28 (2002) 2505 

[25] A.Blondel and F. Jaquet, Proc. ep Facility for Europe, Amsterdam 1979, DESY 79/48 (1979) 391 
(unpublished) 

[26] J.Bliimlein and M.Klein, Snowmass Workshop 1990, proceedings, ed. R.Craven, (1990) 549 
[27] U. Bassler and G. Bernardi, Nucl Instr. Meth. A 361 (1995) 197 

[28] S. Bentvelsen et al., HERA Workshop 1991, Proceedings, eds. W. Buchmiiller and G.Ingelman, 
VI (1992) 23 

[29] K.Hoeger, HERA Workshop 1991, Proceedings, eds. W. Buchmiiller and G.Ingelman, VI (1992) 

43 

[30] ZEUS Collaboration, M. Derrick et al, Z. Phys. C 72 (1996) 399 

[31] A. Glazov, DIS05 Workshop, April 2005, Madison, WI, USA, proceedings, eds. W.Smith and 
S.R. Dasu, (2005) 237 

[32] HI Collaboration, I. Abt et al., Nucl Instr. and Meth. A 386 (1997) 310 Ibid. 386 (1997) 348 



67 



[33] ZEUS Collaboration, "The ZEUS Detector, Status Report 1993", DESY 1993 (unpublished) 

[34] I. Antoniadis, L. Baulieu and C. Kounnas, Nucl. Phys. B 168 (1980) 394 . 

[35] M. Klein and T. Riemann, Z. Phys. C 24 (1984) 151 

[36] E. Derman, Phys. Rev. D 7 (1973) 2755 

[37] C. Callan and D. Gross, PRL 22 (1969) 156 

[38] C. Y. Prescott at al., Phys. Lett. B 84 (1979) 524 

[39] S. L. Glashow, Nucl. Phys. 22 (1961) 579 

A. Salam, Proceedings of the 8th Nobel Colloquium, ed. Svartholm, N. (Almqvist and Wiksell, 
Stockholm 1968) 

S. Weinberg, PRL 19 (1967) 1264 

[40] A. Zee, F. Wilczek and S.B. Treiman, Phys. Rev. D 10 (1974) 2881 
G. AltareUi and G. MartineUi, Phys. Lett. B 76 (1978) 89 

[41] A. Arbuzov, D. Y. Bardin, J. Bluemlein, L. Kalinovskaya and T. Riemann, Comput. Phys. Com- 
mun. 94 (1996) 128 

[42] Particle Data Group, W.M. Yao et al, J. Phys. G 33 (2006) 1 

[43] Inclusive cross section papers of the ZEUS Collaboration: 
M. Derrick et al., Phys. Lett. B 303 (1993) 183 
M. Derrick et al., Phys. Lett. B 316 (1993) 412 
M. Derrick et al., Z. Phys. C 65 (1995) 379 
M. Derrick et al., Phys. Rev. Lett. 75 (1995) 1006 
M. Derrick et al., Z. Phys. C 69 (1996) 607 
M. Derrick et al., Z. Phys. C 72 (1996) 399 
J.Breitweg et al. Phys. Lett. B 407 (1997) 432 
J.Breitweg et al., Eur. Phys. J. C 7 (1999) 609 
J. Breitweg et al. Phys. Lett. B 487 (2000) 53 

J.Breitweg et al. Eur. Phys. J. C 12 (2000) 411 and Erratum Eur. Phys. J. C 27 (2003) 305 

S. Chckanov ct al, Eur. Phys. J. C 21 (2001) 443 

S. Chekanov et al, Phys. Lett. B 539 (2002) 197 and Erratum Eur. Phys. J. C 552 (2003) 308 
S. Chekanov et al, Eur. Phys. J. C 32 (2003) 1 
S. Chekanov et al, Phys. Rev. D 70 (2004) 052001 
S. Chekanov et al, Phys. Lett. B 637 (2006) 210 
S. Chekanov et al, Eur. Phys. J.C 49 (2007) 523 

[44] Inclusive cross section papers of the HI Collaboration: 
T.Ahmed et al.,Phys. Lett. B 299 (1993) 385 
I.Abt et al., Nucl. Phys. B 407 (1993) 515 
I.Abt et al., Phys. Lett. B 321 (1994) 161 
T.Ahmed et al.,Phys. Lett. B 324 (1994) 241 
T.Ahmed et al., Nucl. Phys. B 439 (1995) 471 
S. Aid et al., Nucl. Phys. B 470 (1996) 3 
S. Aid et al., Phys. Lett. B 379 (1996) 319 
C.Adloff et al., Nucl. Phys. B 497 (1997) 3 
C.Adloffet al., Eur. Phys. J. C 13 (2000) 609 
C.Adloff et al., Eur. Phys. J. C 19 (2001) 269 
C.Adloff et al., Eur. Phys. J.C 21 (2001) 33 
C.Adloff et al, Eur. Phys. J. C 30 (2003) 1 
A. Aktas et al., Phys. Lett. B 634 (2006) 173 



68 



HI Collaboration, "Measurement of the Inclusive ep Scattering Cross Section at low and x at 
HERA", to be submitted for publication, Hamburg 2008. 

[45] A. De Rujula, S. L. Glashow, H. D. Politzer, S. B. Treiman, F. Wilczek and A. Zee, 
Phys. Rev. D 10 (1974) 1649 

[46] A. Donnachie and P.V. Landshoff, Proceedings HERA Workshop, Hamburg 1987, ed. R. Peccei, 
Vol.1, p.351 and CERN-Th.5020/88(1988) 

[47] M.Gliick, E.Reya and A.Vogt, Z. Phys. C 53 (1992) 127 

[48] Y. Watanabe et al., PRL 35 (1975) 898 

[49] K.Prytz, Phys. Lett. B 311 (1993) 286 

[50] A.Mueller, Nucl. Phys.: Proc. Suppl. A 29 (1992) 275 

[51] A.DeRoeck, M.Klein and T.Naumann, Phys. Lett. B 385 (1996) 411 

[52] G. AltareUi, R. D. Ball and S. Forte, Nucl. Phys. B 742 (2006) 1 and references therein. 

[53] D. Haidt, Eur. Phys. J. C 35 (2004) 519 

[54] E. A. Kuraev, L. N. Lipatov and V. S. Fadin, Sov. Phys. JETP 44 (1976) 443 and Sov. Phys. 

JETP 45 (1977) 199 

Y. Y. Balitsky and L. N. Lipatov, Sov. Journ. Nucl. Phys. 28 (1978) 822 
[55] ZEUS Collaboration, J.Breitweg et al., Eur. Phys. J. C 7 (1999) 609 
[56] HI Collaboration, C. Adloff et al, Phys. Lett. B 520 (2001) 183 
[57] A. Donnachie and P.V. Landshoff, Z. Phys. C 61 (1994) 139 
[58] ZEUS Collaboration, J. Breitweg et al. Phys. Lett. B 487 (2000) 53 
[59] P. Newman, Int. J. Mod. Phys. A 19 (2004) 1061 
[60] C. Ewerz and O. Nachtmann, Phys. Lett. B 648 (2007) 279 

[61] C. E. Hyde-Wright, L. Frankfurt, M. Strikman and C. Weiss, DIS06 Workshop, April 2006, 
Tsukuba, Japan, proceedings, eds. M. Kuze, K. Nagano and K. Tokushuku, (2006) 279 

[62] P. Jacobs and X.-N. Wang, Prog. Part. Nucl. Phys. 54 (2004) 443 

[63] S. Kretzer, F.I. Olness, R.J. Scalise, R.S. Thorne and U.K. Yang, Phys. Rev. D 64 (2001) 033003 
M. Nakahata, to appear in the proceedings of LP07, Aug. 2007, Daegu, Korea 

[64] L.A.T. Bauerdick, A. Glazov and M. Klein, HERA Physics Workshop, 1996/97, proceedings, eds. 
G. Ingelman, A. DcRocck and R. Klanner VI (1996) 77 

M. Klein, DIS04, Strbske Pleto, 1994, proceedings, eds. D. Bruncko, J. Ferencei and P. Strizenec 
VI (1994) 309 

J. Feltesse Ringberg Workshop, October 2005, proceedings http://www.mppmu.mpg.de/ 

[65] S. Moch, J.A.M. Vermaseren and A. Vogt, Phys. Lett. B 606 (2005) 123 ; R. Thorne, Ringberg 
Workshop, October 2005, proceedings http://www.mppmu.mpg.de/ 

[66] HI Collaboration, C. Adloff et al, Phys. Lett. B 393 (1997) 452 ; see also Eur. Phys. J. C 13 (2000) 
609 

[67] HI Collaboration, C. Adloff et al., Eur. Phys. J.C 21 (2001) 33 . 

[68] HI Collaboration, Paper submitted to EPS2003, July 2003, Aachen, Germany HI prel. 03-043, and 
T. Lastovicka, EPJC direct, Eur. Phy. J. C 33 (2004) s01,s388-s390 

[69] HI Collaboration, "Running at Low Proton Beam Energies" , Expression of Interest, subm. to the 
DESY PRC, October 2005, PRC 11/05, Hl-10/05-622, (unpublished) 

[70] Yu. L. Dokshitzcr, Sov. Phys. JETP 46 (1977) 641 

[71] V. N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438,675 



69 



[72] G. Altarelli and G. Parisi, Nucl. Phys. B 126 (1977) 298 

[73] J. A. M. Vermaseren, A. Vogt and S. Moch, Nucl. Phys. B 691 (129) 2004 and Nucl. Phys. B 724 

(2005) 3 

[74] For detailed discussions and further references see: Small x Collaboration, B. Andersson et al., 
Eur. Phys. J. C 25 (2002) 77 and Eur. Phys. J.C 35 (2004) 67 
J. R. Andersen et al, Eur. Phys. J. C 48 (2006) 53 

[75] M. Ciafaloni, Nucl. Phys. B 296 (1988) 49 

S. Catani, F. Fiorani and G. Marchesini, Phys. Lett. B 234 (1990) 339 
S. Catani, F. Fiorani and G. Marchesini, Nucl. Phys. B 336 (1990) 18 

G. Marchesini, Nucl. Phys. B 445 (1995) 49 

[76] J. B. Dainton, M. Klein, P. Newman, E. Perez and F. Willeke, JINST 1 (2006) PlOOOl ; see also 
http : //www . Ihec . org . uk/ 

[77] M. Botje, Phys. J. C 14 (2000) 285 

http : //www . nikhef . nl/~h24/qcdnuin/ 

[78] C. Pascaud and F. Zomer, preprints LAL-94-42 (1994) and 95-05 (1995) (unpublished) 

[79] MRST, now MSTW Collaboration, see for example: 

A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Eur. Phys. J. C A (1998) 463 
A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Eur. Phys. J. C 23 (2002) 73 
A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Phys. Lett. B 652 (2007) 292 

[80] CTEQ Collaboration, see for example: 

H. L. Lai et al., Eur. Phys. J. C 12 (2000) 375 
J. Pumplin et al., JHEP 0207 (2002) 012 

H. L. Lai, et al. JHEP 0704 (2007) 089 

[81] S. Alekhin, Phys. Rev. D 68 (2003) 014002 

S. Alekhin, K. Melnikov and F. Petriello, Phys. Rev. D 74 (2006) 054033 

[82] E 866/NuSea Collaboration, R. S.Towell et al., Phys. Rev. D 64 (2001) 052002 

[83] M.Klein and B. Reisert, HERA-LHC Workshop, March 2004 to March 2005, proceedings, ed. 
A. DeRoeck and H. Jung, VI (2005) 99 

[84] See the PDF section of the HERA-LHC Workshop, March 2004 to March 2005, proceedings, ed. 
A. DeRoeck and H. Jung, VI (2005) 41 

[85] For a recent paper see: 

W.K. Tung, et al., JHEP 0702 (2007) 053 and references therein 

[86] R. S. Thorne and R. G. Roberts, Phys. Rev. D 57 (1998) 6871 

[87] ZEUS Collaboration, S. Chekanov et al., Z. Phys. C 21 (2001) 441 

[88] ZEUS Collaboration, S. Chekanov et al, Phys. Rev. D 67 (2003) 012007 

[89] HI CoUaboration, C. Adloff et al., Eur. Phys. J.C 30 (2003) 1 

[90] A.D.Martin et al, Eur. Phys. J.C 23 (2002) 73 

[91] J. Pumplin et al., JHEP 0207 (2002) 012 

[92] B. Portheault, PhD Thesis, 03/05, Univ. Paris XI Orsay (2005) (unpublished) 
available at HI http://www-hl.desy.de 

[93] A. Vogt in DIS99, Proc. Suppl. Nucl. Phys. B 79 (1999) 102 , ed. J. Bliimlein and T. Riemann 

[94] ZEUS Collaboration, Phys. Rev. D 67 (2003) 012007 

[95] M. Klein, to appear in the Proceedings of DIS07, April 2007, Munich, Germany 
[96] S. Catani and M.H. Seymour, Nucl. Phys. B 485 (1997) 291 



70 



[97] E. Mirkes and D. Zeppenfeld, Phys. Lett. B 380 (1996) 205 

[98] D. Graudenz, arXiv:hep-ph/970244. 

[99] A.D. Martin et al., Phys. Lett. B 531 (2002) 216 

[100] J. Pumplin et al., JHEP 0207 (2002) 012 

[101] S.V Chekanov, the 14th Topical Conference on Hadron Collider Physics, Sept. 2002, Karlsruhe, 
proceedings, eds. M. Erdmann and T. Mueller (2003) 478. ArXiv:hep-ph/0211298vl. 

[102] G.C. Blazcy ct al.. Physics at RUN II: QCD and Weak Boson Physics Workshop, Nov. 
1999, Batavia, proceedings, eds. U. Baur, R.K. Ellis and D. Zeppenfeld (2000) 47. ArXiv:hep- 
ex/0005012v2. 

[103] S. Catani, et al., Nucl. Phys. B 406 (1993) 187 

[104] S.D. Ellis and D.E. Soper, Phys. Rev. D 48 (1993) 3160 

[105] G. Ingelman, A. Edin and J. Rathsman, Comp. Phys. Comm 101 (1997) 108 

[106] L. Lonnblad, Comp. Phys. Comm 71 (1992) 15 

[107] L. Lonnblad, Z. Phys. C 65 (1995) 285 

[108] B. Andersson, Phys. Rep. 97 (1983) 31 

[109] S. Prixione and G. Ridolfi, Nucl. Phys. B 507 (1997) 315 

[110] G. Corcclla ct al, JHEP 0101 (2001) 010 

G. Marchesini et al., Comp. Phys. Comm 67 (1992) 465 

[111] T. Sjostrand et al., Comp. Phys. Comm 135 (2001) 238 
T. Sjostrand et al, Comp. Phys. Comm 82 (1994) 74 

[112] HI Collaboration, C. Adloff et al., Eur. Phys. J. C 13 (2000) 397 
ZEUS Collaboration, J. Breitweg et al., Phys. Lett. B 479 (2000) 37 
ZEUS Collaboration, S. Chekanov et al., Eur. Phys. J. C 35 (2004) 487 
HI Collaboration, A. Aktas et al., Eur. Phys. J. C 37 (2004) 141 

[113] P. Aurenche, J. Guillet and M. Fontannaz, ZPC 64 (1994) 621 

[114] R. Brun et al.. Technical Report CERN-DD/EE/84-1, CERN, (1987) 

[115] ZEUS Collaboration, S. Chekanov et al., Phys. Lett. B 547 (2002) 164 

[116] HI Collaboration, A. Aktas ct al., DESY-07-073, submitted to Phys. Lett. B 

[117] R.P. Feynman, Photon- Hadron Interactions, Benjamin, New York, 1972 
K.H. Streng, T.F. Walsh and P.M. Zerwas, Z. Phys. C 2 (1979) 237 

[118] A.D. Martin et al., Eur. Phys. J. C 4 (1998) 463 
A.D. Martin et al., Eur. Phys. J. C 14 (2000) 133 

[119] ZEUS Collaboration, S. Chekanov et al., Eur. Phys. J. C 23 (2002) 615 

[120] M. Gliick, E. Reya and A. Vogt, Phys. Rev. D 45 (1992) 3986 

[121] M. Gliick, E. Reya and A. Vogt, Phys. Rev. D 46 (1992) 1973 

[122] HI Collaboration, A. Aktas et al., Phys. Lett. B 639 (2006) 21 . 

[123] ZEUS Collaboration, S. Chekanov et al., DESY-07-092, accepted by Phys. Rev. D 

[124] ZEUS Collaboration, S. Chekanov et al., Nucl. Phys. B 765 (2007) 1 

[125] HI Collaboration, C. Adloff et al., Phys. Lett. B 515 (2001) 17 

[126] ZEUS Collaboration, S. Chekanov et al., Eur. Phys. J. C 44 (2005) 183 

[127] ZEUS CoUaboration, S. Chekanov et al., Phys. Lett. B 649 (2007) 12 



71 



[128] ZEUS Collaboration, S. Chekanov et al., DESY-07-102, accepted by Nucl. Phys. B 
ZEUS Collaboration, S. Chekanov et al., DESY-07-100, accepted by Eur. Phys. J. C 
ZEUS Collaboration, S. Chekanov et al, DESY-07-062, accepted by Nucl. Phys. B 
ZEUS Collaboration, S. Chekanov et al, Phys. Lett. B 649 (2007) 12 
HI Collaboration, A. Aktas et al., Eur. Phys. J. C 51 (2007) 271 
ZEUS Collaboration, S. Chekanov et al., Nucl. Phys. B 729 (2007) 492 
HI CoUaboration, A. Aktas et al., Eur. Phys. J. C 46 (2006) 27 
HI Collaboration, A. Aktas et al., Eur. Phys. J. C 37 (2004) 141 
HI Collaboration, A. Aktas et al., Eur. Phys. J. C 33 (2004) 477 

[129] ZEUS Collaboration, J. Breitweg et al., Phys. Lett. B 507 (2001) 70 
ZEUS CoUaboration, S. Chekanov et al., Phys. Lett. B 547 (2002) 164 
HI Collaboration, C. Adloff et al., Eur. Phys. J. C 19 (2001) 289 
ZEUS Collaboration, S. Chekanov et al., Phys. Rev. D 67 (2003) 012007 
ZEUS Collaboration, S. Chekanov et al, Eur. Phys. J. C 42 (2005) 1 
HI Collaboration, C. Adloff et al., Eur. Phys. J. C 21 (2001) 33 
ZEUS CoUaboration, S. Chekanov et al., Phys. Lett. B 558 (2003) 41 
ZEUS CoUaboration, S. Chekanov et al., Eur. Phys. J. C 31 (2003) 149 
ZEUS CoUaboration, S. Chekanov et al., Phys. Lett. B 560 (2003) 7 
ZEUS CoUaboration, S. Chekanov et al, Eur. Phys. J. C 44 (2005) 183 
ZEUS CoUaboration, S. Chekanov et al, Nucl. Phys. B 700 (2004) 3 
ZEUS CoUaboration, S. Chekanov et al, Phys. Lett. B 207 (2006) 12 

[130] C. Glasman, arXiv:hep-ex/0506035 

[131] S. Bethke, arXiv:hep-ex/0407021 

[132] S. Bethke, arXiv:hep-ex/0606035 

[133] HI and ZEUS CoUaborations, Papers submitted to EPS07, July 2007, Manchester, UK 

[134] ZEUS CoUaboration, S. Chekanov ct al., Eur. Phys. J. C 42 (2005) 1 

[135] A.M. Cooper-Sarkar, Ringbcrg Workshop, October 2005, 
proceedings http : // www . mppmu . mpg . de/ 

[136] B.W. Harris and J. Smith, Phys. Rev. D 57 (1998) 2306 

[137] R.S. Thorne, J. Phys. G 25 (1999) 1307 

[138] R.G. Roberts and R.S. Thorne, Eur. Phys. J. C 19 (2001) 339 

[139] M.A.G. Aivazis et al., Phys. Rev. D 50 (1994) 3102 

M.A.G. Aivazis, F.I. Olness and W-K. Tung Phys. Rev. D 50 (1994) 3085 

[140] B.A. Kniehl et al., Phys. Rev. D 71 (2005) 014018 
B.A. Kniehl et al, Eur. Phys. J. C 41 (2005) 199 

[141] P.M. Nadolsky et al, Phys. Rev. D 67 (2003) 074015 

[142] M.G. Bowler Z. Phys. C 11 (1981) 169 

[143] V.G. Kartivehshvih, A.K. Likhoded and V.A. Petrov Phys. Lett. B 78 (1978) 615 

[144] C. Peterson et al, Phys. Rev. D 27 (1983) 105 

[145] M. Ciafaloni, Nucl. Phys. B 296 (1988) 49 

S. Catani, F. Fiorani and G. Marchesini, Phys. Lett. B 234 (1990) 339 
S. Catani, F. Fiorani and G. Marchesini, Nucl. Phys. B 336 (1990) 18 
G. Marchesini, Nucl. Phys. B 445 (1995) 49 

[146] H. Jung and G.P. Salam, Eur. Phys. J. C 19 (2001) 339 

[147] H. Jung, arXiv:hep-ph/0411287 



72 



[148] W.K. Tung, S. Kretzer and C. Schmidt, J. Phys. G 28 (2001) 983 

[149] HI Collaboration, A. Aktas et al., Eur. Phys. J. C 51 (2007) 271 

[150] B.A. Kniehl ct al., Phys. Rev. Lett. 96 (2006) 012001 

B.A. Kniehl, G. Kramer and M. Maniatis, Nucl. Phys. B 711 (2005) 345 Erratum-ibid. 720 (2005) 
345. 

[151] ZEUS Collaboration, S. Chekanov et al., Phys. Rev. D 69 (2004) 012004 

[152] HI Collaboration, C. Adloff et al., Phys. Lett. B 528 (2002) 199 

[153] HI Collaboration, A. Atkas et al., Eur. Phys. J. C 50 (2007) 251 

[154] S. Prixione, P. Nason and G. Ridolfi, Nucl. Phys. B 454 (1995) 3 

S. Prixione, M.L. Mangano, P. Nason and G. Ridolfi, Phys. Lett. B 348 (1995) 633 

[155] G. Heinrich and B.A. Kniehl, Phys. Rev. D 70 (2004) 094035 

[156] ZEUS Collaboration, S. Chekanov et al., Nucl. Phys. B 729 (2005) 492 

[157] ZEUS Collaboration, S. Chekanov et al., Phys. Rev. D 70 (2004) 012008 

[158] HI Collaboration, A. Aktas et al., Eur. Phys. J. C 41 (2005) 453 

[159] ZEUS Collaboration, S. Chekanov et al., Phys. Lett. B 599 (2005) 173 

[160] HI Collaboration, A. Aktas et al., Phys. Lett. B 621 (2005) 56 

[161] ZEUS Collaboration, S. Chekanov et al., Eur. Phys. J. C 50 (2007) 299 

[162] ZEUS Collaboration, Papers submitted to the EPS07, July 2007, Manchester, UK 

[163] A.V. Lipatov and N.P. Zotov, Phys. Rev. D 75 (2007) 014028 

[164] D. Pitzl et al., Nucl. Instr. Meth. A 454 (2000) 334 

[165] E.N. Koffeman, Nucl. Instr. Meth. A 473 (2001) 26 

[166] HI Collaboration, A. Aktas et al., Eur. Phys. J. C 40 (2005) 349 

[167] HI Collaboration, A. Aktas et al., Eur. Phys. J. C 45 (2006) 23 

[168] ZEUS Collaboration, S. Chekanov et al., Phys. Lett. B 599 (2004) 173 

[169] ZEUS Collaboration, S. Chekanov et al., JHEP 07 (2007) 074 

ZEUS Collaboration, S. Chekanov et al., Eur. Phys. J. C 44 (2005) 351 
HI Collaboration, A. Aktas et eil.,Eur. Phys. J. C 38 (2005) 447 
ZEUS Collaboration, J. Breitweg et al., Phys. Lett. B 481 (2000) 213 

[170] ZEUS Collaboration, S. Chekanov et a\.,Phys. Lett. B 649 (2007) 111 

[171] HI Collaboration, A. Aktas et al.,Eur. Phys. J. C 47 (2006) 597 

[172] ZEUS Collaboration, M. Derrick et al., Phys. Lett. B 346 (1995) 399 

[173] HI Collaboration, T.Ahmed et al. Phys. Lett. B 348 (1995) 681 . 

[174] J. Colhns, Phys. Rev. D 57 (1998) 3051 Erratum-ibid. 61 (2000) 019902 

[175] J. Bliimlein and D. Robaschik, Phys. Rev. D 65 (2002) 096002 

[176] For a review see J. R. Forshaw, R. Sandapen and G. Shaw, JHEP 0611 (2006) 025 

[177] HI Collaboration, A. Aktas et al., Eur. Phys. J.C 48 (2006) 715 Ibid 48 (2006) 749 

[178] ZEUS Collaboration, S. Chekanov et al.. Euro. Phys. J. C 52 (2007) 813 

[179] ZEUS Collaboration, S. Chekanov et al., Nucl. Phys. B 713 (2005) 3 

[180] J.Lukasi, to appear in the Proceedings of DIS07, April 2007, Munich, Germany 

[181] HI Collaboration, A. Aktas et al., JHEP 0710 (2007) 042 

[182] A.Martin, M. Ryskin and G.Watt, Phys. Lett. B 644 (2007) 131 

[183] HI Collaboration, A. Aktas et al., Eur. Phys. J. C 51 (2007) 549 . 



73 



[184] The LHC FP420 Project, see http://www.fp420.com/ 

[185] ZEUS Collaboration, S. Chekanov et al., Nucl. Phys. B 713 (2005) 3 

[186] J. Bartel et al., Eur. Phys. J. C 7 (1999) 443 

[187] B.K. Kopeliovitch et al., Phys. Lett. B 324 (1994) 469 

[188] S.J. Brodsky et al., Phys. Rev. D 50 (1994) 3134 

[189] L.L. Frankfurt, A. Freund and M. Strikman, Phys. Rev. D 58 (1998) 114001 Erratum-ibid 59 

(1999) 119901 

D. Miiller et al., Fortschr. Phys. 42 (1994) 347 
X. Ji, Phys. Rev. Lett. 78 (1997) 610 

K. Golec-Biernat and A. Martin, Phys. Rev. D 59 (1999) 014029 

[190] ZEUS Collaboration, S. Chekanov et al., Nucl. Phys. B 695 (2004) 3 

[191] ZEUS Collaboration M. Derrick et al., Eur. Phys. J. C 6 (1999) 603 
HI Collaboration, A. Aktas et al., Eur. Phys. J. C 13 (2000) 371 
HI Collaboration, C. Adloff et al., Phys. Lett. B 483 (2000) 360 
ZEUS Collaboration, J. Breitweg et al., Phys. Lett. B 487 (2000) 273 
ZEUS Collaboration, S. Chekanov et al, Nucl. Phys. B 718 (2005) 3 
HI Collaboration, A. Aktas et al, Eur. Phys. J. C 46 (2006) 585 
ZEUS Collaboration, S. Chekanov et al, PMC Physics A 1 (2007) 6 

[192] ZEUS Collaboration, S. Chekanov et al., Phys. Lett. B 573 (2003) 46 

[193] HI Collaboration, A. Aktas et al, Eur. Phys. J. C 44 (2005) 1 

[194] A. Freund and M.F. McDermott, Phys. Rev. D 65 (2002) 091901 
A. Freund and M.F. McDermott, Eur. Phys. J. C 23 (2002) 651 

[195] A. Freund, M.F. McDermott and M. Strikman, Phys. Rev. D 67 (2003) 036001 

[196] C. Marquet and L. Schoeffel, Phys. Lett. B 639 (2006) 471 

[197] K. Golec-Biernat and M. Wiisthoff, Phys. Rev. D 59 (1999) 014017 

[198] H. Spiesberger et al., in Physics at HERA, Proceedings for the Workshop, Hamburg, 1991, V2 
(1991) 195 

[199] ZEUS Collaboration, S. Chekanov et al., Eur. Phys. J. C 32 (2003) 1 
HI Collaboration, C. Adloff et al., Eur. Phys. J. C 30 (2003) 1 

ZEUS Collaboration, S. Chekanov et al, Phys. Rev. Lett. 539 (2002) 197 ibid. 552 (2003) 308 
HI Collaboration, C. Adloff et al, Eur. Phys. J. C 19 (2001) 269 

[200] ZEUS Collaboration, S. Chekanov et al., Phys. Rev. D 70 (2004) 052001 
ZEUS Collaboration, S. Chekanov et al, Eur. Phys. J. C 28 (2003) 2 

[201] V. Chekehan, to appear in the Proceedings of ICHEP06, July 2006, Moscow. 

[202] HI Collaboration, A. Aktas et al, Phys. Lett. B 634 (2006) 173 

[203] CDF and DO collaborations, ArXiv:hep-ex/0703034. 

[204] J.Bluemlein, M.Klein, T. Naumann and T. Riemann, Workshop on HERA Physics, 1986/7, pro- 
ceedings, cd. R. Peccei, VI (1987) 67; 

R. Beyer et al.. Future HERA Physics, Proceedings for the Workshop, Hamburg, 1995/6, eds. 
R.Klanner and A.De Roeck, VI (1996) 140. 

[205] Y.D. Ri, to appear in the Proceedings of ICHEP06, July 2006, Moscow, Russia 

[206] E. Rizvi and T. Sloan, Eur. Phys. J.direct C 3 (2001) N2 

[207] M. Kuze and Y. Sirois, Prog. Part. Nucl. Phys 50 (2003) 1 Erratum-ibid. 53 (2004) 583 



74 



[208] W. Buchmiiller, R. Riickl and D. Wyler, Phys. Lett. B 191 (1987) 442 Erratum-ibid. 448 (1999) 
320 

[209] J. Butterworth and H. Dreiner, Nucl. Phys. B 397 (1993) 3 

[210] ZEUS Collaboration, S. Chckanov et al., Phys. Rev. D 68 (2003) 052004 

[211] HI Collaboration A. Aktas et al., Phys. Lett. B 629 (2005) 9 

[212] L3 Collaboration, M. Acciarri et al., Phys. Lett. B 486 (2000) 81 

[213] OPAL Collaboration G. Abbiendi et al.,Eur. Phys. J. C 6 (1999) 1 

[214] DO Collaboration, V.M. Abazov et al., Phys. Rev. D 71 (2005) 071104 

[215] DO Collaboration, B. Abbot et al, Phys. Rev. Lett. 79 (1997) 4321 
DO Collaboration, B. Abbot et al., Phys. Rev. Lett. 80 (1998) 2051 
DO Collaboration, V.M. Abazov et al., Phys. Rev. D 64 (2001) 029004 
DO Collaboration, V.M. Abazov et al., Phys. Rev. Lett. 88 (2002) 191801 

[216] ZEUS Collaboration, S. Chckanov ct al., Eur. Phys. J. C 44 (2005) 463 

[217] HI Collaboration A. Aktas et al., DESY-07-009, Submitted to Eur. Phys. J. C 

[218] HI Collaboration, V. Andreev et al., Phys. Lett. B 561 (2003) 241 

[219] ZEUS Collaboration, S. Chekanov et al., Phys. Lett. B 559 (2003) 153 

[220] Y. De Boer, to appear in the Proceedings of DIS07, April 2007, Munich, Germany 

K. Korcsak-Gorzo, to appear in the Proceeding of DIS07, April 2007, Munich, Germany 

[221] HI and ZEUS Collaborations, Paper submitted to LP07, August 2007, Daegu, Korea, Hlprelim- 
07-162, ZEUS-prel-07-029 

[222] ZEUS Collaboration, S. Chekanov et al., Phys. Lett. B 583 (2004) 41 

[223] HI Collaboration, A. Aktas et al, Eur. Phys. J. C 48 (2006) 699 

[224] HI Collaboration, A. Aktas et al, Eur. Phys. J. C 33 (2003) 9 

[225] HI Collaboration, A. Aktas et al, Eur. Phys. J. C 31 (2003) 17 

[226] O. Ota, to appear in the Proceedings of DIS07, April 2007, Munich, Germany 
G. Brandt, to appear in the Proceedings of DIS07, April 2007, Munich, Germany 

[227] HI and ZEUS Collaborations, Paper submitted to EPS07, Manchester, UK, Hlprehm-07-166, 
ZEUS-prel-07-024. 



75 



